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PREFACE 


This  research  was  performed  to  provide  some  quantitative 
estimates  of  the  Army's  attack  helicopter  vulnerability  to 
a hostile  artillery  threat.  To  the  extent  specified  by 
the  stated  assumptions  and  constraints,  the  answers  to  this 
vulnerability  issue  are  available  from  the  simulation  model 
described  in  this  thesis. 

Thanks  are  due  to  Major  Ed  Duff  of  the  Fighter  Division, 
Air  Force  Studies  and  Analyses,  for  his  continuing  interest 
in  this  thesis  and  for  his  many  useful  suggestions  which 
properly  defined  the  scope  of  this  research.  Major 
Dick  Kulp  was  extremely  helpful  in  all  matters  pertaining 
to  probability  theory,  and  the  difficulties  encountered  in 
writing  this  model's  target  location  error  algorithm  were 
resolved  by  him.  Captain  Jerry  Sullivan  was  instrumental 
in  gathering  the  much  needed  feedback  from  the  Army  pilots 
at  Ft.  Hood,  Texas. 

A very  special  thanks  goes  to  my  wife,  Sandy,  for  her 
encouragement  throughout  the  life  of  this  project. 

Emil  H.  Koenig,  III 


(This  thesis  typed  by  Sharon  A.  Gabriel) 
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uniformly  distributed  area  fire  and  precision  fire  against 
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the  modeled  activities  from  one  replication  to  another. 
Current  concepts  of  attack  helicopter  employment  in  the 
anti-armor  role  are  used.  Concepts  of  intervisibility, 
detectability,  target  location  error,  C-cubed  time  delay, 
and  artillery  round-to-round  dispersion  are  also  incorporated 
in  the  model.  The  model  was  experimented  upon,  using  a 
hypothetical  artillery  weapon  system,  and  the  vulnerability 
predictions  of  the  model  are  displayed  and  analyzed. 
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A SIMULATION  MODEL 


OF  ATTACK  HELICOPTER  VULNERABILITY 
TO  HOSTILE  ARTILLERY  FIRE 

I . INTRODUCTION 

This  thesis  is  a modeling  effort  to  predict  the  vulner- 
ability of  the  Army's  attack  helicopter  (AH)  to  hostile 
artillery  fire  while  the  helicopter  is  being  operated  near 
the  Forward  Edge  of  the  Battle  Area  (FEBA)  in  its  armor 
defeating  role.  The  research  reported  here  involves  the 
development  of  a computer  simulation  model  that  calculates, 
within  appropriate  assumptions  and  constraints,  the  AH 
probability  of  kill  (P^)  in  this  particular  battlefield 
environment . 

To  simulate  the  armor  defeating  role,  the  AH  movement 
on  the  battlefield  is  structured  to  conform  to  the  Army's 
current  concept  of  employment  for  this  weapon  system.  The 
effect  o'  battlefield  threat  systems  other  than  hostile 
artillery,  whether  acting  independently  or  cumulatively,  is 
not  simulated.  Additionally,  "near  the  FEBA"  is  taken  to 
mean  a discrete  portion  of  the  land  mass  within  which 
artillery  is  the  only  simulated  threat  to  AH  survival,  and 
outside  of  which  no  threat  to  the  AH  exists.  The  hostile 
artillery  fire  is  laid  down  in  two  formats:  (i)  a uniform 
area  fire  simulating  the  enemy's  preparation  of  an  area  of 
intended  conquest,  and  (ii)  precise  indirect  fire  simulating 


observer-directed  action  against  an  AH  point  target  of 
high  interest. 

The  predicted  vulnerability  during  these  specified 
battlefield  conditions  is  a primary  accomplishment  of  this 
work.  The  sensitivity  of  this  vulnerability  to  changes  in 
several  of  the  main  variables  of  interest  is  also  of 
importance  because  of  the  amount  of  uncertainty  present  in 
estimating  parameter  values. 

The  Problem  Statement 

The  problem  is  that  very  little  quantitative  information 
exists  which  predicts  AH  vulnerability  to  hostile  artillery 
fire  while  the  helicopter  is  being  operated  near  the  FEBA. 

The  shallow  extent  to  which  this  topic  has  been 
previously  studied  was  verified  by  a rather  thorough 
literature  search  of  documents  held  by  both  the  Defense 
Documentation  Center  and  the  Defense  Logistics  Studies 
Information  Exchange.  The  vulnerability  mentioned  in  this 
problem  statement  is  defined  to  mean  the  probability  that 
the  AH,  during  any  single  exposure  event  in  the  target  area, 
will  experience  an  aircraft  kill  from  any  or  all  of  the 
effects  of  the  nearby  detonation  of  high  explosive,  contact 
detonating  artillery  warheads.  As  used  in  this  thesis,  an 
aircraft  kill  is  either  attrition  of  the  helicopter  or  a 
forced  landing  within  30  minutes . 
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Thesis  Arrangement 

Because  some  classified  data  were  used  in  this  study, 
a thesis  format  was  necessary  which  would  allow  for  the 
widest  dissemination  of  the  methodology  and,  yet,  safeguard 
the  sensitive  information.  As  a result,  this  thesis  is 
written  in  two  documents,  the  bulk  of  which  is  this  main 
report  that  includes,  by  way  of  example,  the  effects  of  a 
hypothetical  artillery  system.  When  specific  enemy  threat 
data  versus  the  AH  were  considered,  the  results  are  discussed 
in  the  classified  Appendix  B to  this  main  report.  When  it 
becomes  necessary  to  acknowledge  classified  documents  in 
this  main  report,  the  reference  will  be  by  title  only. 

Goals 

Several  goals  were  set  at  the  outset  of  this  work  in 
order  to  properly  measure  the  degree  to  which  the  thesis' 
overall  objective  was  being  met,  that  objective  being  to 
develop  a method  to  answer  "how  vulnerable  is  the  AH  to 
hostile  artillery  fire?".  These  goals  were: 

1.  Construct  a computer  model  to  simulate  the 
artillery  fire  versus  AH  environment. 

2.  Test  that  this  model  behaves  properly  by  applying 
both  internal  verification  and  external  validation. 

3.  Perform  sensitivity  analysis  upon  the  main 
parameters  of  interest. 

4.  Draw  general  conclusions  about  the  AH  vulnerability 
to  the  artillery  threat  based  upon  the  model's  output. 
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Importance  of  This  Study  ! 

There  are  several  reasons  why  this  problem  is  being 
studied  in  addition  to  the  purely  academic  goal  of  gaining  i 

a better  understanding  of  what  computer  simulation  is  all  \ 

about.  Foremost,  this  is  a topic  of  current  Air  Force 
interest.  The  Fighter  Division,  Directorate  of  Studies 

I 

and  Analyses,  Headquarters  Air  Force,  suggested  this  topic 

/ 

as  a candidate  for  a thesis  study  and  this  office  has 
actively  supported  the  research.  Secondly,  by  constructing 
a computer  model  and  then  documenting  the  thought  processes 
used  in  the  model's  development,  one  should  obtain  a better 
understanding  of  the  real  world  physical  processes  whereby 
artillery  is  indeed  a threat  of  some  concern  to  the  AH 

force.  Finally,  by  experimenting  with  the  model  over  some 

r 

} reasonable  ranges  of  the  interesting  parameters,  an  under- 

standing may  be  gained  as  to  how  sensitive  helicopter 

I 

vulnerability,  as  predicted  by  this  model,  is  to  changes 
in  friendly  tactics  and  enemy  capabilities. 

Model  Application 

If  this  simulation  model  succeeds  in  addressing  the 
relevant  issues  and  interactions  among  battlefield  * 

resources,  it  may  find  broad  use  in  developing  input 
vulnerability  data  to  tactical  simulation  models  of  far 
larger  scope.  One  example  of  a larger  scope  model  is  the 
USAF  TAC  WARRIOR  model.  It  is  also  clear  that  there  is 
yet  much  to  be  accomplished  in  studying  this  relatively 
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small  facet  of  the  overall  Close  Air  Support  issue. 
Hopefully,  this  thesis  can  act  as  the  basis  model  for 
any  future  sophistication  that  may  be  desired. 

Thesis  Overview 

Chapter  II,  Background,  puts  this  thesis  in  proper 
perspective  by  reviewing  previous  studies  which  have 
addressed  this  problem.  Included  in  this  chapter  is  a 
description  of  the  various  weapon  systems  that  were  modeled 
in  this  simulation.  Chapter  III,  The  Methodology,  explains 
what  real  world  events  were  made  a part  of  the  model  and  the 
manner  in  which  these  events  were  allowed  to  interact.  Also 
in  Chapter  III  is  the  discussion  of  how  the  model  was 
validated.  Chapter  IV,  The  Simulation  Model,  contains  the 
details  cf  the  computer  model  including  a description  of 
the  model's  logic  and  mechanics.  Chapter  V,  Sample  Results, 
demonstrates  the  model's  capability  by  showing  the  output 
produced  from  the  consideration  of  a hypothetical  artillery 
weapon.  Finally,  Chapter  VI,  Summary,  Conclusions,  and 
Recommendations , provides  a recap  of  the  main  points,  draws 
conclusions  based  upon  the  model's  performance,  and  illumi- 
nates those  shortcomings  in  the  present  effort  which,  if 
corrected,  should  produce  a better  product. 
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II.  BACKGROUND 


This  chapter  is  intended  to  provide  the  reader  with  a 
better  perspective  of  the  scope  of  this  problem.  It  does 
this  by  describing  some  important  points  made  by  previous 
studies  of  this  topic  and  related  issues.  Following  this 
are  some  observations  which  have  been  made  about  the  AH's 
vulnerability  to  artillery  fire.  Finally,  the  primary 
weapon  systems  which  are  modeled  in  this  computer  simulation 
are  described  in  sufficient  detail  to  acquaint  the  reader 
with  their  capabilities  and  concepts  of  employment. 

Previous  Studies 

Many  models  have  been  constructed  which  attempted  to 
quantify  the  events  which  happen  in  a land  battle  and  the 
subsequent  effects  of  those  events  upon  the  forces  in  the 
field.  A large  share  of  this  effort  are  the  models  that 
attempt  to  simulate  the  exchange  of  all  types  of  fires 
between  opposing  forces  and  to  specify  the  result  of  such 
fire  on  the  disposition,  strength  and  combat  capability, 
and  the  viable  options  which  remain  for  the  engaged  forces. 
One  subset  of  these  models  includes  fairly  detailed  studies 
of  the  Army's  AH  as  it  is  being  used  in  the  role  for  which 
it  is  now  designed  : anti-armor.  It  certainly  appears 
that  such  in  depth  analysis  is  fully  warranted  because  it 
has  been  noted  by  Taylor  that  "the  U.  S.  Army  has  to  hold 
the  Fulda  Gap  with  fewer  than  fifty  Huey  Cobra  armed 
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helicopters,  against  thousands  of  Warsaw  Pact  tanks...." 

(Ref.  1:29). 

STATE  II  is  a computer  simulation  which  models  the 
armor  versus  anti-armor  battle  at  the  company  level  (Ref.  2). 
It  was  used  as  the  starting  point  for  this  thesis  because 
it  addressed  a firepower  issue  and  it  structured  the  real 
world  problem  into  a computer  program  that  was  highly  event 
oriented.  Although  helicopters  were  not  specifically 
mentioned  as  being  players  in  this  STATE  II  model,  the  AH 
concept  and  its  inherent  anti-armor  capability  could  be 
adapted  to  this  simulation. 

Another  report  presented  the  results  of  an  Army  field 
experiment  in  which  the  primary  purpose  was  to  investigate 
"scout  helicopter  effectiveness  in  one  of  a variety  of 
missions  on  the  mid- intensity  battlefield"  (Ref.  3:1-1-1). 

An  interesting  observation  made  in  this  report  was  that, 

* 

if  the  helicopter  exposed  itself  during  an  unmask  maneuver 
for  a period  of  37  seconds,  this  "represented"  a 0.15 
probability  that  a ground  delivered  projectile  would  impact 
near  the  scout  vehicle.  This  probability  relationship  was 
not  substantiated  with  an  explanation  of  the  type  of 
projectile  (although  direct-fire  tank  weapon  could  be 
inferred)  or  of  how  this  relationship  was  derived. 

* 

The  unmask  maneuver  is  a flight  maneuver  whereby  the 
helicopter  leaves  its  hidden  position  behind  cover  or 
concealment  in  a pop-up  to  an  altitude  sufficient  to  provide 
line  of  sight  to  the  target  of  interest.  The  remask  is  the 
reciprocal  maneuver. 
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An  older  study  that  described  an  artillery  assessment 
model  which  was  a tool  used  in  a war  game  called  SYNTAC 
considered  artillery  fire  against  helicopters,  but  only 
against  a fixed  percentage  of  the  helicopter  force  that 
would  be  in  the  artillery's  target  area.  The  difficulty 
with  this  approach  is  that  this  fixed  percentage  was  the 
approximation  of  the  numbers  of  helicopters  that  would  be 
considered  parked  or  on  ground  alert  (Ref.  4:13).  By 
design,  this  simulation  model  would  not  address  the  issue 
of  artillery  against  a moving  AH.  It  appears  that  this 
study  assumed  that  the  artillery  is  no  threat  at  all  to 
a helicopter  in  motion. 

There  have  also  been  several  recent  classified  Army 
studies  which  examined  in  great  detail  the  artillery  threat 
to  the  AH.  These  models  differ  from  the  model  developed  in 
this  thesis  in  the  level  of  detail  of  the  environment  which 
is  incorporated  into  the  simulation.  For  instance,  one 
particular  Army  model  requires  that  the  user  input  a detailed 
description  of  the  terrain  over  which  the  AH  will  fly,  thus 
the  results  of  such  a model  are  applicable  only  to  that 
particular  geographic  region.  This  thesis,  by  assuming 
average  terrain  features,  presents  a model  that  is  more 
universally  applicable.  Also,  this  same  Army  model  inputs 
specific  point  and  area  targets  for  the  artillery  fire 
missions  to  engage.  These  inputs  were  made  as  a result  of 
close  scrutiny  of  the  terrain  map  by  artillery  officers. 

By  contrast,  this  thesis  simulates  area  wide  coverage  of 
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the  artillery  fire  and,  in  addition,  simulates  the  forward 
observer  (FO)  who  attempts  to  direct  fire  specifically 
against  the  AH  when  detection  occurs. 

The  Army  is  generally  concerned  about  the  AH  vulner- 
ability to  artillery  fire.  As  was  discovered  through 
dialogue  with  personnel  at  Ft.  Hood,  Texas,  the  artillery 
threat  is  viable  enough  to  have  influenced  the  planned 
disposition  of  the  AH  forces  in  the  target  area  (Ref.  5). 

This  concern  is  taken  into  account  by  the  present  tactics 
of  AH  employment.  More  simply  stated  is  the  fact  that 
"quite  obviously  the  helicopter  has  to  vacate  the  artillery 
barrage  areas,  as  will  all.  other  thin-skinned  vehicles" 

(Ref.  6:57). 

The  AH-1Q/S 

This  Army  attack  helicopter  is  being  produced  in  two 
versions,  models  Q and  S,  both  of  which  are  derivatives  of 
the  original  Army  Huey  Cobra  AH-1G.  The  main  difference 
between  these  newer  versions  and  the  basic  aircraft  is  that 
the  newer  versions  have  been  equipped  to  employ  the  TOW 
(Tube-launched,  Optically-tracked,  Wire-guided)  anti-tank 
missiles.  This  potent  addition  to  the  Army's  anti-armor 
arsenal  is  generally  referred  to  as  the  Cobra/TOW. 

This  AH  is  a single  engine  aircraft  driving  a two  bladed 
rotor.  It  is  operated  by  a crew  of  two  whose  stations  are 
arranged  in  tandem.  The  gunner  occupies  the  forward 
station,  the  pilot  occupies  the  rearward  station.  This 
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weapon  system  is  reported  to  have  a maximum  level  speed  of 
123  knots  with  missile  launchers  on  board  and  a maximum 
range  of  507  km  (Ref.  7:274-275). 
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The  primary  armament  carried  by  the  AH-1Q/S  is  the 
TOW  missile.  This  missile  carries  a high  explosive  shaped 
charge  warhead  which  was  specifically  designed  to  penetrate 
armor.  The  gunner  provides  guidance  commands  to  the  weapon 
from  the  time  of  launch  until  target  impact  occurs  by  keeping 
the  target  image  centered  in  his  gyro-stabilized  gun  sight. 
These  guidance  signals  are  passed  to  the  missile  through  two 
wires  which  are  unreeled  during  the  missile's  entire  fly  out. 
One  AH-1Q/S  can  carry  up  to  eight  missiles  at  once.  The 
maximum  range  attributed  to  this  TOW  missile  system  is 
3750  meters  (Ref.  8). 

The  role  of  the  AH  in  modern  battle,  as  stated  in 
Field  Manual  17-50,  is  to  "destroy  sufficient  numbers  of 
the  enemy  to  convince  him  to  break  off  his  attack,  to  give 
up  a defensive  area,  or  to  move  away  from  an  area  vital  to 
friendly  forces"  (Ref.  9:4-1).  The  method  of  destruction  is, 
of  course,  the  employment  of  the  TOW  missiles.  However, 
the  AH  does  not  intend  to  accomplish  this  mission  by  acting 
alone.  His  partners  on  the  battlefield  are  the  Aeroscout 
helicopter  and  the  friendly  ground  forces,  both  of  which 
find  tanks  and  other  armored  vehicles  and  provide  this 
information  to  the  AH.  The  Aeroscout  helicopter  is  a more 
capable  teammate  for  the  AH  than  surface  forces  because  the 
helicopter's  high  degree  of  mobility  insures  more  frequent 
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contact  with  targets  of  the  opposing  side.  In  addition, 
the  Aeroscout  assumes  the  responsibility  of  providing  local 
protection  to  the  AH  while  he's  engaging  the  targets. 

A sequence  of  events  that  the  AH  could  be  expected  to 
perform  as  he  works  with  the  Aeroscout  on  a combat  mission 
might  be: 

Move  to  holding  area 

Coordinate  with  Aeroscout 

Move  to  battle  position  and  receive 
target  hand  off 

Partially  unmask 

Acquire  target 

Unmask  as  required  to  fire 

Engage  - Remask  - Move  to  alternate 
firing  positions  or  return  to  holding 
area  (Ref.  9:43) 

A holding  area  is  a site  that  the  AH  occupies,  usually  only 
for  brief  periods,  while  the  Aeroscout  is  busy  coordinating 
other  AH  movements  into  attack  positions. 

Threat  Artillery 

As  is  well  known,  the  purpose  of  employing  artillery 
is  primarily  to  reinforce  the  offensive  striking  power  of 
an  aggressor.  This  is  typically  characterized  in  Soviet 
doctrine  by  the  employment  of  massive  fires,  if  the 
artillery  resources  are  available  (Ref.  10:2-13).  The 
big  bore  weapons  usually  remain  under  the  division  or 
regiment  operational  control,  but  are  supportive  of  the 
fire  and  maneuver  plans  of  subordinate  units  as  necessary. 
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As  an  example  of  this  Combined  Arms  Armies  concept  of 
warfare,  it  is  quite  possible  that  a motorized  rifle  or 
tank  battalion  will  receive  the  support  of  six  large 
artillery  pieces  during  an  offensive  operation  (Ref.  10:2-19). 

Two  common  artillery  threat  systems  are  the  122  mm 
Towed  Howitzer  (D-30)  and  the  152  mm  Towed  Gun-Howitzer 
(D-20).  The  122  mm  weapon  is  reported  to  have  a maximum 
range  of  15,300  meters  and  a rate  of  fire  of  seven  to  eight 
rounds  per  minute.  The  152  mm  weapon  is  credited  with  a 
17,000  meter  range  and  four  rounds  per  minute  rate  of  fire 
(Ref.  11:2-19).  The  probability  of  kill  data  for  these 
two  threat  systems  against  the  basic  model  AH-1G  and  the 
result  of  the  simulation  process  while  incorporating  these 
specific  weapon  characteristics  are  contained  in  the 
classified  Appendix  B. 

Forward  Observer 

The  forward  observer  (FO)  , as  simulated  in  this  model, 
is  a combatant  who  is  equipped  with  an  armored  vehicle 
which  permits  the  FO  to  remain  near  the  other  combat 
vehicles  at  the  FEBA.  His  primary  function  is  to  see  as 
much  of  the  battlefield  situation  as  is  possible  and  to 
report  sightings  of  targets  of  interest  such  as  the  Army's 
AH.  The  goal  of  this  FO  is  to  detect  the  presence  of  an 
AH  anytime  that  intervisibility  (an  unobstructed  line  of 
sight)  exists  between  his  position  and  the  AH's  position, 
to  measure  as  accurately  as  is  practicable  the  location  of 
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the  AH,  and  to  radio  this  target  sighting  to  an  artillery 
unit  so  that  a concentrated  artillery  fire  mission  may  be 
brought  to  bear  upon  that  sighting.  Specific  Soviet  FO 
capability  and  equipment  is  not  included  in  this  model. 
Instead,  variables  of  interest  which  impact  upon  the  FO's 
performance  are  programmed  within  a wide  range  of  values  to 
permit  parametric  study. 


Ill . THE  METHODOLOGY 
Simulation  Requirements 

There  are  only  a few  basic  requirements  that  the 
simulation  model  must  perform  in  order  that  it  become 
representative  of  the  real  world  environment.  The  most  > 

important  requirement  is  the  proper  synthesizing  of  the 
AH's  movement  on  the  battlefield.  The  desired  movement 
algorithm  must  be  one  that  is  conceptually  correct  for  the 
environment  being  modeled.  For  this  thesis,  the  AH  was 
allowed  enough  freedom  of  play  to  insure  adequate  randomness 
from  one  replication  to  the  next;  yet,  this  randomness 
was  sufficiently  bounded  so  that  the  overall  pattern  of 
movement  is  representative  of  the  Army's  current  employment 
doctrine . 

Another  important  modeling  requirement  is  that  the 
artillery  be  required  to  fall  in  the  same  area  in  which  the 
AH  is  operating.  Randomness  is  incorporated  into  this 
function  both  in  the  uniform  area  fire  and  in  the  precision 
fire  against  a point  target.  « 

The  last  major  requirement  is  that  simulated  time  j 

within  the  model  must  be  permitted  to  remain  "frozen" 
while  the  effects  of  a critical  event,  such  as  an  artillery 
impact  near  the  AH  location,  is  properly  assessed. 
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C3  delays,  and  weapon  inaccuracies.  The  fact  that  such 
artillery  support  is  readily  available  is  consistent  with 
the  idea  that  a relatively  stable  situation  is  being 
modeled  as  opposed  to  such  fast  moving  battlefield  operations 
as  exploitation. 

SIMSCRIPT  11.5 

The  SIMSCRIPT  II. 5 computer  language  was  selected  as 
the  language  of  this  computer  model  for  several  reasons. 

The  main  reason  is  that  SIMSCRIPT  II. 5 is  appropriately 
structured  to  handle  time  dependent  routines  or  events. 

The  language's  built-in  timing  routine  greatly  simplified 
the  scheduling  of  the  many  events  which  the  model  is 
performing.  Another  useful  feature  of  this  language  is 
that  an  interface  is  provided  which  allows  for  linking  a 
FORTRAN  subroutine,  as  necessary,  to  the  main  simulation 
program.  This  particular  ability  was  required  to  enable 
the  execution  of  the  plotting  package  which  produced  many 
of  the  figures  contained  in  this  thesis. 

Because  of  software  incompatibility,  the  AFIT 
implementation  of  the  SIMSCRIPT  II. 5 language  did  not 
allow  for  the  user  of  the  model  to  be  interactive  with 
this  program  at  an  on-line  terminal.  This  mild  limitation 
presented  no  difficulty  in  the  development  of  this  thesis; 
however,  interactive  capability  may  be  worthwhile  in  future 
attempts  to  update  or  change  this  model. 
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Goals  Accomplished 


The  first  goal  of  this  thesis  was  accomplished  by  the 
construction  of  the  computer  model  described  in  this  paper. 
This  was  done  by  building  into  this  model  the  known 
characteristics  of  the  weapon  systems,  a typical  scenario, 
and  an  established  employment  doctrine. 

The  second  goal  was  the  testing  of  the  behavior  of 
the  model.  This  was  conducted  in  two  ways.  The  internal 
verification  was  done  to  insure  that  the  functions  within 
the  model  were  operating  as  planned.  This  proper  behavior 
was  confirmed  by  the  examination  of  numerous  plots  of  the 
model's  data  output.  These  plots  were  executed  on  AFIT's 
CDC  3292  Incremental  Plotter  device.  An  example  of  the  type 
of  graphics  which  were  produced  by  the  plotting  package  is 
contained  in  Figure  1. 

This  is  a typical  ground  track  portrayal  of  the  AH 
movement  which  begins  at  A,  proceeds  to  firing  positions 
close  to  the  "front"  of  the  operating  area,  then  exits 
at  B on  the  return  leg  to  re-arm  and/or  re-fuel.  The 
plotting  routine  which  accomplished  this  is  called 
SUBROUTINE  SCRIBE  and  its  FORTRAN  code  listing  is  included 
in  Appendix  A along  with  the  SIMSCRIPT  II. 5 code  for 
this  computer  model. 

The  external  validation  of  the  model  was  done  by 
preparing  a series  of  representative  graphics,  similar  to 
Figure  1,  which  were  then  sent  to  Army  personnel  for  informal 
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Figure  1.  Sample  Graphics 

critique.  The  personnel  who  participated  in  this  feedback 
of  information  were  highly  qualified  Cobra/TOW  aviators, 
usually  with  Vietnam  combat  experience,  assigned  to  the 
6th  Air  Cavalry  Combat  Brigade  at  Ft.  Hood,  Texas.  An 
interesting  example  of  the  kinds  of  information  that  were 
provided  by  the  Army  included  a strong  recommendation  that 
the  AH  holding  area  not  be  included  within  the  initial 
target  area  boundaries  of  3000  meters  by  3000  meters.  This 
holding  area  was  originally  built  into  the  model  to  conform 
to  Army  guidance  in  Field  Manual  17-50  which  envisioned 
the  holding  area  as  being  a location  near  enough  to  the  FEBA 
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to  allow  for  the  AH's  fast  response  to  any  target  assignment. 
However,  as  pointed  out  by  the  Cobra/TOW  pilots  at  Ft.  Hood, 
this  holding  area  would  be  selected  based  upon  the  reality 
of  the  situation  and  this  includes  avoidance  of  the  enemy's 
artillery  fires,  if  possible.  As  a result,  the  holding 
area  concept  was  taken  out  of  the  model. 

The  model  was  also  validated  to  a certain  degree  by 
the  adaptation  of  two  algorithms  that  were  obtained  from 
outside  sources.  One  algorithm  was  the  line  of  sight  model 
which  was  written  by  David  C.  Hardison,  currently  an 
Under  Secretary  of  the  Army,  and  included  in  a paper  presented 
to  the  Terrain  Modeling  Working  Group  at  a Colorado  Springs 
meeting  in  August  1977.  It  has  been  used  in  the  joint 
Air  Force/Army  project  called  the  NATO  Fire  Support  Require- 
ments Study . 

The  other  algorithm  is  the  detection  subroutine  of 
the  Night  Combat  Model  developed  by  the  Army's  Night  Vision 
Laboratory.  This  program  has  been  used  as  the  acquisition 
probability  subroutine  for  the  Army  Materiel  Systems 
Analysis  Activity  model  called  War  Game  and  the  Battalion 
Level  Differential  Model.  The  Army's  Concepts  and  Analysis 
Agency  is  also  using  this  algorithm  in  the  Carmonette  Model. 
Both  routines  are  fully  explained  in  Chapter  IV. 

The  third  goal,  which  was  the  sensitivity  analysis 
upon  the  main  parameters  of  interest,  was  performed  by 
making  many  single-factor  experiments  of  the  model.  The 
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impact  upon  the  of  varying  the  values  of  the  main 
parameters  are  shown  in  graphical  form  in  Chapter  V for  the 
hypothetical  artillery  system  and  in  Appendix  B for  the 
Soviet  122  mm  and  152  mm  Howitzers. 

The  final  goal  was  to  make  general  conclusions  about 
AH  vulnerability  to  the  artillery  threat.  These  conclusions 
were  made  based  upon  the  output  from  this  computer  model. 
The  conclusions  are  a part  of  Chapters  V and  VI. 

Primary  Variables  of  Interest 

The  user  of  this  model  is  required  to  select  the 
values  of  many  different  variables;  thus,  any  particular 
run  or  series  of  runs  will  be  expressly  tailored  to  the 
user's  needs.  Each  variable  value,  of  course,  can  become 
the  basis  for  as  much  in  depth  testing  of  the  model's 
results  as  is  desired.  Through  discussion  with  personnel 
at  the  Fighter  Division,  Directorate  of  Studies  and  Analyses, 
Headquarters  Air  Force,  it  was  decided  that  four  of  the 
variables  in  this  model  would  be  examined  in  some  detail. 

One  variable  is  the  amount  of  time  that  the  AH  remains 
in  its  initial  battle  position  prior  to  commencing  the  first 
TOW  engagement.  This  was  built  into  the  model  to  simulate 
the  time  required  for  final  strike  coordination  between 
the  AH  and  the  Aeroscout.  In  this  model,  that  amount  of 
time  is  represented  by  the  input  variable  CALL. TO. FIRE. 

This  loiter  time  is  not  applied  to  any  subsequent  reattacks 
that  the  AH  may  accomplish  during  any  replication.  This 
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simulates  that  additional  targets  are  seen  by  the  AH 
during  its  initial  TOW  engagement  and  that  the  Aeroscout 
is  not  needed  to  generate  these  subsequent  targets.  This 
CALL. TO. FIRE  variable  was  allowed  to  take  on  values  between 
0 seconds  and  180  seconds.  It  should  be  noted  that  a fixed 
time  delay  is  incorporated  at  each  attacking  position, 
whether  it's  the  initial  battle  position  or  a reattacking 
position.  This  time  constant  is  35  seconds,  which  is  the 
written  guidance  to  the  AH  crew,  and  it  represents  the  upper 
bound  of  time  which  should  be  spent  in  the  TOW  firing 
operation  (Ref.  9:4-9). 

Another  variable  is  the  intensity  of  the  artillery's 
area  fire.  Since  the  area  fire  continues  without  slack  or 
interruption  during  the  AH's  entire  stay  in  the  target 
area,  the  intensity  is  dependent  only  upon  the  rate  of  fire. 
Control  of  this  intensity  parameter  is  achieved  by  simply 
changing  the  value  of  T1  which  is  the  time  interval,  in 
seconds,  between  successive  artillery  round  impacts.  This 
T1  parameter  was  allowed  to  have  values  between  2 seconds 
and  20  seconds. 

The  target  location  error  was  another  variable  of 
interest.  This  parameter  is  the  radius,  in  meters,  of  a 
circle  drawn  about  the  AH's  exact  position.  Inside  of 
this  circle,  all  points  are  equally  likely  to  be  selected 
as  the  FO's  determination  of  where  the  AH  truly  lies. 

Labeled  TLE,  this  target  location  error  takes  on  values 
between  10  meters  and  100  meters. 
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The  last  variable  of  interest  in  this  thesis  is  the 
time  which  elapses  between  the  instant  that  the  FO  detects 
the  presence  of  the  AH  target  and  the  instant  that  the 
first  round  of  the  precision  fire  hits  the  target  area. 

The  concept  of  this  time  delay  embodies  such  mission  related 
activities  as  : the  FO  must  translate  this  target  detection 
into  an  appropriate  set  of  map  coordinates,  the  request 
for  the  fire  mission  must  be  radioed  to  the  supporting 
artillery  unit,  the  artillery  unit  must  convert  the  target 
coordinates  to  the  appropriate  tube  alignment  parameters, 
and  the  first  round  must  experience  the  appropriate  time 
of  flight.  It  was  desired  that  this  timing  delay  variable 
would  primarily  represent  the  C3  delay  aspect  of  the  above 
chain  of  activities;  thus,  it  is  named  C3. DELAY  and  it 
ranged  in  value  between  0.5  minutes  and  3 minutes. 

Replication  Planning  Factors 

To  compile  the  object  deck  on  AFIT's  CDC  6600  system 
takes,  on  the  average,  19  seconds  of  Central  Processing  (CP) 
time  and  requires  71000  octal  words  of  Central  Memory. 

Trial  runs  have  shown  that  the  execution  phase  of  the  model 
generally  consumes  14.4  seconds  of  CP  time  per  100  replica- 
tions. One  replication  is  one  entry  and  exit  of  the  AH 
for  the  Area  of  Operations,  including  all  AH  movement  in 
between,  all  artillery  activity,  and  any  damage  assessment 
calculations  which  may  be  necessary. 
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It  was  desired  that  the  overall  P^  calculation  be 
within  ±0.03  of  the  true,  but  unknown,  Pj, . Also,  a 
90%  level  of  confidence  was  desired  for  the  range  selected. 
As  shown  in  Shannon  (Ref.  12:192),  the  number  of  required 
replications  (n)  to  achieve  the  0.03  upper  and  lower 
P^  bound  (d)  is  related  to  the  two-tailed  standardized 
normal  statistic  (Z)  by 


n 


(1) 


where  Z ,,,  = 1.645  for  the  90%  confidence  level, 

a / 2 

This  yields  752  required  replications,  which  equates 
to  approximately  108  seconds  of  CP  execution  time  for  each 
P^  data  point  desired.  Based  upon  this  analysis,  the  model 
was  replicated  752  times  for  each  changing  level  of  a main 
parameter.  This  plan  was  followed  for  all  of  the  computed 
data  points  which  are  plotted  in  Chapter  V and  in  Appendix  B. 


23 


IV.  THE  SIMULATION  MODEL 


This  chapter  explains  the  development  of  the  individual 
routines  of  the  simulation  model.  The  routines  of  the 
SIMSCRIPT  II. 5 language  are  labeled  as  events  and  these 
events  are  scheduled  for  execution  by  the  built-in  timing 
routine.  The  complete  code  listing  for  this  model  is 
contained  in  Appendix  A. 

The  AH  is  always  constrained  to  remain  within  the  Area 
of  Operations  (AO)  boundary,  whose  size  is  determined  by 
the  user.  These  dimensions,  labeled  X.MAX  and  Y.MAX,  were 
set  at  values  of  3000  meters  each  for  all  of  the  sample 
replications  in  this  research.  The  AO  is  depicted  in 
Figure  2. 

The  x-axis  of  this  coordinate  system  is  the  line 
farthest  away  from  the  armor  threat , and  the  y = 3000  meters 
line  is  the  closest.  Conceptually,  the  enemy's  first  line 
of  armored  vehicles  is  simulated  to  be  on  a line  at 
y = 5500  meters.  Thus,  it  is  assumed  that  the  AH  does  not 
proceed  any  closer  than  2500  meters  to  the  armored  threat 
which  the  AH  is  designed  to  defeat. 

Event  NOE 

Event  NOE,  when  scheduled,  generates  a randomly 
oriented  flight  path  between  the  AH's  currently  occupied 
position  and  the  ultimate  destination.  This  is  accomplished 
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Figure  2.  Area  of  Operations 

by  generating,  at  fixed  time  intervals,  individual  velocity 
vectors  for  the  AH  to  maintain.  These  small  velocity 
vectors  randomly  vary  from  one  to  another  in  magnitude  and 
direction.  This  was  done  to  simulate  nap-of -the-earth 
flight  which  Field  Manual  17-50  defines  as  flying  a 
"weaving  slalom-like  route  within  the  planned  corridor 
while  remaining  oriented  along  the  general  axis  of  movement" 
(Ref.  9:4-8).  These  velocity  vectors  were  used  to  specify 
the  series  of  discrete  points  that  the  AH  is  to  occupy  on 
its  travels  in  this  particular  mode  of  flight.  The  flight 
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path  maintained  between  each  of  these  discrete  points  was 
assumed  to  be  straight  line  and  at  constant  velocity. 

The  magnitude  of  each  velocity  vector  was  a function 
of  the  fixed  time  interval  (DEL.T)  and  a variable  airspeed 
term.  This  variable  airspeed  was  randomly  drawn  from  a 
uniform  distribution  with  an  upper  and  lower  bound  of 
LOSPEED  and  (0 . 5 ) LOSPEED , respectively.  The  user  inputs 
the  DEL.T  value  in  seconds  and  the  LOSPEED  value  in  knots. 

For  this  thesis,  values  of  6 seconds  and  15  knots  were 
used . 

The  direction  of  the  vectors  also  varied  from  one  to 
another.  This  part  of  the  algorithm  consisted  of  defining 
the  straight  line  path  between  the  AH's  currently  occupied 
position  and  the  ultimate  destination.  An  angular  deviation 
from  this  path  was  randomly  drawn  from  a uniform  distribution 
between  ± PHI. MAX.  PHI. MAX  is  entered  by  the  user  in  units 
of  degrees,  and  45  degrees  was  the  value  used  in  this 
research . 

Since  the  AH's  currently  occupied  position  is  always 
stored  by  the  model  and  the  magnitude  and  direction  of  the 
newly  generated  velocity  vector  has  been  specified,  vector 
addition  was  used  to  provide  the  coordinates  of  the  AH’s 
next  position.  This  iterative  process  was  executed  once 
each  6-second  period  until  the  AH  finally  arrived  at  its 
destination.  An  example  of  the  nap-of -the-earth  ground  track 
generated  by  the  Event  NOE  is  shown  in  Figure  3. 
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RUN  *6 


Figure  3.  Nap-of-the-Earth  Ground  Track 

All  of  the  movement  which  is  depicted  above  the 
y = 2000  meters  line  was  generated  by  the  Event  NOE.  The 
constraint  that  the  AH  must  employ  nap-of -the-earth  flight 
whenever  it  is  within  1000  meters  of  the  front  edge  of 
the  AO  was  formulated  through  discussion  with  the  Ft.  Hood 
aviators  previously  mentioned  in  this  thesis. 

Event  LOS 

Event  LOS  determines  whether  or  not  unobstructed 
line  of  sight  (intervisibility)  exists  between  the  FO  and 
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the  AH.  Event  LOS  is  represented  by  the  top  step  of  the 
logic  flow  chart  in  Figure  4.  An  adaptation  of  the 
Hardison  line  of  sight  algorithm,  first  mentioned  in 
Chapter  III,  is  used.  This  event  is  excited  once  every 
10  seconds  of  simulated  time.  The  calculation  of  the 
probability  that  line  of  sight  exists  is  a function  of 
target -to-observer  range  ( R ) , terrain  type  coefficient 

( K ) , and  the  observer's  viewing  height  ( H ) above 

the  terrain.  These  variables  are  related  to  the  probability 
of  line  of  sight  by: 


f Li  • 

where  is  K | A^  + (H/100)  1 J , R is  expressed 

in  kilometers,  H is  in  meters,  and  A.  , B.  , C. 
are  modified  target  coefficients.  For  this  thesis,  values 
of  At  = 1.56  , Bi  = 3.48  , Ci  = 1.173  , K = 1.0  , 
and  H = 3 meters  were  used.  The  A^  , B^  , and 
coefficients  were  sample  values  suggested  by  the  Fighter 
Division,  Directorate  of  Studies  and  Analyses,  Headquarters 
Air  Force.  The  K term  is  based  upon  the  type  of  terrain 
being  simulated  and  it  has  these  representative  values: 


smooth  terrain 

K = 

1.9 

rolling  terrain 

K = 

1.0 

rough  terrain 

K = 

0.3 
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(2) 


The  user  of  the  model  must  specify  the  values  of  K and  H 


This  Event  LOS,  after  calculation  of  the  P^Qg  at 
each  10-second  interval,  draws  a random  deviate  from  a 
uniform  ( 0,1  ) distribution  and  compares  this  deviate 
value  with  P^Qg.  For  those  instances  when  the  deviate  value 
is  less  than  or  equal  to  P^os*  ■'■^ne  sight  is  presumed  to 
exist  and  the  detection  algorithm  (Event  DETECT)  is  then 
executed.  Otherwise,  another  line  of  sight  iteration  is 
scheduled  after  the  required  10-second  delay. 


Event  DETECT 

Event  DETECT  is  essentially  a language  modified  edition 
of  the  Army's  Night  Vision  Laboratory-developed  Detection 
Subroutine  of  the  Night  Combat  Model.  The  original  FORTRAN 
coding  was  rewritten  in  SIMSCRIPT  II. 5 language.  This 
algorithm  is  executed  for  the  target -to -observer  range  at 
the  instant  that  line  of  sight  exists  in  the  model.  The 
algorithm  calculates  the  probability  of  target  detection 
as  a complex  function  of  the  viewing  device's  capability 
to  respond  to  the  target's  presence  and  of  the  observer's 
capability  to  recognize  this  device  response  as  a target 
signature.  The  routine  calculates  this  detection  probability 
by  applying  light  levels  and  target/background  contrasts 
to  the  number  of  resolvable  line  scan  cycles  required  by  the 
observer  to  make  a target  detection.  A curve-fitting 
technique  is  then  employed  to  fit  empirical  data  to  these 
specific  encounter  parameters  to  arrive  at  the  appropriate 
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probability  of  detection. 

The  following  assumptions  were  built  into  the  Event 
DETECT.  The  FO  must  acquire  and  be  able  to  classify  a 
target  image  to  a sufficient  degree  as  to  be  able  to 
distinguish  between  tracked  versus  wheeled  vehicle.  The 
meteorological  conditions  are  seven  kilometer  visibility 
range  and  the  atmospheric  attenuation  coefficient  is  0.558  . 
The  reflective  contrast  between  AH  and  its  background 
conditions  is  a fixed  ratio  of  0.3  . The  sky  brightness 

to  ground  brightness  ratio  is  3.0  . The  FO  is  searching 

a 20°  x 30°  area  and  using  7-power  binoculars.  Thirty 
percent  of  the  available  light  into  the  binoculars  is  lost 
due  to  optical  attenuation.  The  only  user  input  into  this 
algorithm  is  the  ambient  light  level  (AL1),  measured  in 
foot -candles  , at  the  target.  This  research  used  AL1  = 1000 
foot-candles  to  simulate  bright  sunlight. 

The  end  product  of  this  algorithm  is  a computation  of 
the  probability  that  the  FO  will  detect  the  target  within 
a 10-second  period  after  line  of  sight  conditions  occur. 

This  value  is  called  P10.  A random  deviate  is  then  drawn 
from  a uniform  ( 0,1  ) distribution  and  this  value  is  compared 
to  P10  to  make  a decision  as  to  whether  or  not  detection  has 
occurred.  If  detection  occurred,  a precision  fire  artillery 
mission  is  begun;  otherwise,  the  line  of  sight  event  is 
re-scheduled . 
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Event  PREC .FIRE 


Event  PREC. FIRE  determines  the  magnitude  of  the 
precision  fire  activity.  The  user  controls  this  process  by 
specifying  the  maximum  number  of  rounds  that  may  be  employed 
during  any  single  precision  fire  mission.  This  parameter 
is  called  ROUNDS.  Also,  the  user  specifies  the  value  of 
ROF , which  is  the  rate  of  fire  parameter.  Values  of 
ROUNDS  = 40  and  ROF  = 8 rounds  per  minute  were  used  in 
this  research. 

While  a precision  fire  mission  is  being  executed,  the 
line  of  sight  event  in  the  model  is  discontinued.  Event 
PREC. FIRE  re-initiates  the  line  of  sight  algorithm  after 
the  precision  fire  mission  has  been  completed. 

Event  PREC. FIRE  also  incorporates  the  concept  of  a 
C3  time  delay  between  when  the  AH  detection  has  occurred 
and  when  the  first  round  of  the  volley  impacts  the  target 
area.  The  user  specifies  this  value  of  C3. DELAY  in 
minutes  and,  for  this  thesis,  a value  of  one  minute  was 
used  for  those  replications  which  did  not  apply  any  sensi- 
tivity to  this  variable. 

The  algorithm  which  derived  the  weapon's  aiming 
point  coordinates  is  also  contained  in  Event  PREC. FIRE. 

At  the  instant  of  AH  detection,  the  AH's  position  (X22,Y22) 
is  passed  to  this  event.  It  was  desired  that  the  forward 
observer  (FO)  specify  an  aiming  point  which  would  randomly 
lie  within  a circle  of  radius  TLE,  centered  about  (X22,Y22). 


Also,  this  randomness  concept  was  interpreted  to  mean  that 
all  points  within  the  "TLE"  circle  could  be  selected  by 
the  FO  as  the  ultimate  aiming  point  with  equal  probability. 

To  achieve  this  condition,  a polar  coordinate  frame  of 
reference  was  used  with  the  origin  at  (X22,Y22)  and  the 
initial  line  positioned  as  would  a positively  oriented 
x-axis.  The  direction  coordinate  (BETA)  was  randomly 
drawn  from  a uniform  distribution  between  0 and  2 n radians. 
The  length  coordinate  (RHO)  is  the  product  of  the  user's 
specified  target  location  error  (TLE)  and  the  square  root 
of  a random  deviate  drawn  from  a uniform  (0,1)  distribution. 
This  BETA  and  RHO  deviation  was  then  applied  to  the  (X22,Y22) 
coordinate  to  specify  the  aiming  coordinate  (AIMX,AIMY). 

For  those  replications  of  the  model  which  did  not  vary  the 
target  location  error,  a fixed  value  of  TLE  = 50  meters 
was  used. 

Event  THREAT 

Event  THREAT  serves  many  functions.  If  the  model  has 
scheduled  this  event  in  order  to  apply  another  artillery 
impact  for  the  simulated  area  or  barrage  fire  mission,  the 
required  impact  coordinate  is  drawn  from  two  different 
uniform  distributions  whose  range  of  values  are  from  0.0 
to  X.MAX  and  from  0.0  to  Y.MAX.  This  was  done  to  make  any 
point  within  the  AO  as  equally  likely  to  be  selected  as 
any  other  point.  The  user  controls  the  intensity  of  this 
artillery  activity  by  specifying  the  value  of  T1 . This  T1 
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variable  is  the  time  interval,  in  seconds,  between  successive 
artillery  impacts.  Figure  5 is  the  plotting  of  all  area 
fire  impacts  that  occurred  during  a single  replication  of 
the  model  in  which  the  AH  remained  in  the  AO  for  8.409 
minutes.  The  frequency  of  artillery  impacts  was  once  every 
five  seconds. 


Figure  5.  Area  Fire 

However,  if  the  Event  THREAT  is  scheduled  because 
detection  of  the  AH  has  initiated  a precision  fire  artillery 
mission,  this  routine  determines  the  impact  point  of  the 


round  based  upon  the  aiming  coordinates  (AIMX,AIMY)  which 
were  passed  to  it  from  the  Event  PREC.FIRE.  The  user  has 
control  of  the  footprint  size  of  the  precision  fire  mission 
by  specifying  values  for  PE.R  and  PE.D,  in  meters.  These 
parameters  are  the  probable  error  in  range  and  probable 
error  in  deflection  of  the  single  gun  being  simulated. 

Since  probable  error  is  generally  defined  to  mean  0.6745a 
(a  = one  standard  deviation  of  the  distribution  of  the 
fall  of  shot),  the  algorithm  which  specified  the  random 
deviate  that  represents  the  gun's  dispersion  drew  from  a 
normal  probability  distribution  function  of  mean  0.0  and 
standard  deviation  PE/0.6745.  The  deflection  deviate  was 
added  to  the  x-coordinate  of  the  aiming  point  and  the  range 
deviate  was  added  to  the  y-coordinate . An  example  of  an 
application  of  40  rounds  of  precision  fire  against  an  aiming 
point  is  depicted  in  Figure  6. 

As  is  shown  in  Figure  6,  the  uniformly  laid  area  fire 
continues  to  be  applied,  unaffected  by  the  generation  of 
the  precision  fire  mission.  In  Figure  6,  a rectangle  was 
drawn  which  was  centered  about  the  gun's  aiming  coordinates. 
Based  upon  the  definition  that  PE  = 0.6745a,  on  the  average 
such  a rectangle  should  contain  50%  of  the  impacts  of  this 
precision  fire  algorithm. 

The  final  function  of  Event  THREAT  is  to  compute  the 
P^  of  each  artillery  impact  against  the  AH  and  store  this 
value  for  future  computation.  This  operation  is  accomplished 
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Figure  6.  Precision  Fire  with  Area  Fire 


by  first  checking  the  range  between  AH  position  and 
artillery  impact  at  the  instant  of  simulated  detonation. 

For  those  rounds  which  land  outside  of  the  maximum  range 
at  which  any  hazard  to  the  AH  is  expected,  the  computation 
is  disregarded  and  the  Event  THREAT  is  terminated.  For 
those  rounds  which  land  inside  of  this  maximum  range,  a 
search  is  made  through  the  threat  P^  data  array  and  the  two 
P^  data  entries  which  are  nearest  to  the  range  in  question 
are  found.  A linear  interpolation  between  these  end  values 
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of  is  performed.  The  resultant  for  this  particular 
encounter  is  then  stored  in  an  array  called  PKS  and  the 
artillery's  impact  coordinates  are  stored  in  arrays  called 
PK.X  and  PK.Y.  Tests  made  upon  this  P^  algorithm  have 
shown  that  this  linear  interpolation  technique  never  exceeds 
a 10%  error  from  the  actual  P^  value  and  it  generally  is 
within  3%. 

The  Scenario  Events 

There  are  several  events  in  this  model  which  are 
related  to  one  another  in  the  sense  that  these  events  control 
the  details  of  the  scenario  which  make  any  replication 
different  from  any  other.  Their  event  names  are 
SCENARIO,  INGRESS,  ENTRY. LEG,  1. REATTACK,  2. REATTACK, 

3. REATTACK,  DROPBACK,  and  EGRESS. 

SCENARIO  acts  as  the  traffic  controller  by  scheduling 
other  scenario  related  events,  as  appropriate.  Events 
ENTRY. LEG  and  EGRESS  control  the  two  high  speed  profiles, 
which  simulate  low  level  flight,  in  which  the  AH  is  in  a 
straight  line  flight  path  either  into  or  out  of  the  AO. 

The  speed  at  which  this  maneuver  is  accomplished  is  controlled 
by  the  user  when  he  specifies  the  value  of  the  variable 
HISPEED,  in  knots.  Once  specified  at  the  beginning  of  program 
execution  , HISPEED  remains  constant.  HISPEED  = 50  knots 
was  the  value  used  in  this  research. 

Event  DROPBACK  controls  the  transition  between  the 
nap-of-the-earth  profiles,  which  are  used  by  the  AH  near 
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the  simulated  TOW  firing  positions,  and  the  high  speed 
egress  from  the  AO.  As  an  additional  function,  Event 
ENTRY. LEG  provides  similar  transition  between  the  high 
speed  ingress  and  subsequent  nap-of -the-earth  flight.  The 
coordinates  at  which  these  transitions  are  made  are  labeled 
in  the  model  as  (XSLOWDOWN.YSLOWDOWN)  and  (XSPEEDUP , YSPEEDUP) . 
Events  1. REATTACK,  2. REATTACK,  and  3. REATTACK  control  the 
execution  of  the  AH's  movement  to  alternate  firing  positions, 
as  dictated  by  Event  INGRESS.  It  should  be  noted  that  the 
unmask  and  remask  flight  maneuvers,  usually  associated  with 
TOW  engagements  at  the  firing  positions,  is  not  simulated 
in  this  model. 

Event  INGRESS  is  the  workhorse  of  this  entire  package 
of  events.  The  most  important  product  of  this  event  is  the 
random  determination  of  how  many  reattacks  are  to  be  made 
by  the  AH,  if  any  at  all.  During  any  one  replication,  it 
was  desired  that  four  options  be  available  to  the  AH 
corresponding  to  0,  1,  2,  or  3 reattacks,  and  that  each 
option  be  equally  likely  to  be  selected.  This  was  done 
by  drawing  a random  integer  (12)  from  a uniform  (0,3) 
distribution.  If  12  was  returned  with  a value  other  than 
zero,  any  subsequent  reattacking  locations  were  biased 
toward  the  center  of  the  AO  to  prevent  the  AH  from  ever 
exceeding  the  lateral  boundaries.  Also,  these  reattacking 
locations  were  allowed  to  vary  from  the  previously  occupied 
firing  position  uniformly  between  ± 50  meters  in  the 


y-coordinate  and  uniformly  between  30  and  150  meters  in 
the  x-coordinate . Additionally,  the  Event  INGRESS  specifies, 
in  a random  manner,  the  entry  and  exit  points,  the  initial 
battle  position,  the  FO ' s position,  and  the  SLOWDOWN  and 
SPEEDUP  coordinates . 

Figure  7 is  an  example  of  what  the  AH's  ground  track,  as 
generated  by  this  model,  might  look  like. 


Figure  7.  Sample  Ground  Track,  No  Reattacks 

During  this  particular  replication,  the  AH  was  not 
directed  to  perform  any  reattacks.  Consequently,  after  the 


allotted  time  was  spent  at  the  initial  battle  position, 
the  AH  employed  nap-of -the-earth  flight  until  a point  was 
reached  that  was  1000  meters  away  from  the  forward  edge  of 
the  AO.  Then,  low  level  flight  at  50  knots  was  flown  from 
that  point  along  the  shortest  straight  line  segment  to  the 
rearward  edge  of  the  AO.  The  AH  spent  7.675  minutes 
executing  the  ground  track  which  is  depicted  in  Figure  7.  / 

This  may  be  compared  with  Figure  8,  which  shows  a different 
replication  of  AH  ground  track  activity. 


RUN  # 7 


Figure  8.  Sample  Ground  Track,  3 Reattacks 
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As  can  be  seen,  the  randomly  generated  scenario  within 
the  model  produces  completely  different  ground  tracks  at 
each  replication.  For  the  run  depicted  in  Figure  8,  the 
AH  made  three  reattacks.  The  initial  battle  position  is 
the  upper  right  hand  point  of  the  ground  track.  Because  of 
the  biasing  feature,  subsequent  battle  positions  were 
generated  progressively  closer  to  the  centerline  of  the  AO. 
During  the  replication  shown  in  Figure  8,  the  AH  was 
operating  in  the  AO  for  11.126  minutes. 

Event  SIMEND 

Event  SIMEND  is  the  logical  ending  of  any  replication. 

Because  of  this,  the  event  handles  most  of  the  statistical 

bookkeeping,  including  the  calculation  of  P^.  The  approach 

used  was  that  of  assigning  a of  zero  for  any  replication 

which  generated  no  impacts  close  enough  to  warrant  any 

probability  of  kill.  However,  for  those  instances  in  which 

there  were  at  least  one  impact  encounter  which  warranted 

P^  consideration,  each  impact  was  examined  in  turn  against 

a random  deviate  drawn  from  a uniform  (0,1)  distribution. 

If  a random  deviate  was  less  than  or  equal  to  the  probability 

of  kill  for  any  single  encounter,  an  AH  kill  was  confirmed 

and  the  overall  Pj,  for  that  replication  was  set  at  1.0. 

Otherwise,  the  replication's  P^  was  set  at  0.0.  The  mean 

value  of  the  P,  statistic  was  tallied  under  the  label  STAT.l. 
k 

This  average  P^  per  AH  operation  in  the  AO  is  the  output 
received  by  experimenting  with  this  model. 
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Event  SIMEND  controls  when  the  plotting  package,  called 
SUBROUTINE  SCRIBE,  is  executed.  The  determination  of 
whether  the  number  of  replications,  as  specified  by  the 
user,  has  been  completed  is  also  made  in  this  event.  If 
it  is  determined  that  at  least  one  more  replication  of  the 
model  is  necessary,  global  variables,  which  were  primarily 
used  as  counters,  are  set  to  zero  before  the  next  replication 
is  begun.  Also,  any  events  named  THREAT,  PREC.FIRE,  and/or 
LOS,  which  may  be  waiting  for  execution  in  the  timing  queue, 
are  cancelled  and  destroyed. 

Programs  MAIN  and  SCRIBE 

The  two  main  functions  of  the  program  MAIN  are  to  allow 
the  user  to  input  the  value  of  certain  variables  and  to 
initially  schedule  the  occurrence  of  the  events  called 
SCENARIO,  THREAT,  and  LOS. 

The  subroutine  called  SCRIBE  is  written  in  FORTRAN  code 
and  contains  the  instructions  which  produced  most  of  the 
figures  contained  in  this  thesis.  The  user  of  the  model  is 
cautioned  that  this  particular  FORTRAN  linkage  to  the  basic 
SIMSCRIPT  II. 5 model  allows  for  only  a one-direction  flow 
of  arguments  from  the  calling  routine  (SIMEND)  to  the 
SUBROUTINE  SCRIBE.  However,  the  SIMSCRIPT  II. 5 language 
does  permit  the  passage  of  argument  values  in  both  directions 
if  the  user  so  desires  and  the  appropriate  control  cards 
are  inserted. 
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This  completes  the  detailed  explanation  of  the 
simulation  model.  As  can  be  seen  from  these  descriptions, 
this  event -s tructured  model  is  well  matched  to  the 
capabilities  of  the  SIMSCRIPT  II. 5 language. 


V.  SAMPLE  RESULTS 


To  offer  an  example  of  the  AH  vulnerability  calculations 
that  this  model  provides,  a hypothetical  artillery  threat 
was  created.  The  curve  which  relates  the  versus  range 
to  target  for  this  fictitious  threat  is  graphed  in  Appendix  A. 

The  vulnerability  results  which  are  graphed  in  this 
chapter  are  a function  of  the  changing  levels  of  the  model's 
variables  called  T1  , CALL .TO . F IRE , TLE , and  C3. DELAY.  T1 
is  the  time  interval  between  successive  artillery  impacts  of 
the  simulated  area  fire,  measured  in  seconds.  CALL. TO. FIRE 
is  the  number  of  seconds  that  the  AH  must  remain  in  place, 
after  arrival  at  the  initial  battle  position,  prior  to 
commencing  a TOW  engagement.  TLE  is  the  forward  observer's 
target  location  error,  measured  in  meters.  C3. DELAY  is 
the  amount  of  time  between  the  forward  observer's  detection 
of  the  AH  and  the  first  impact  of  the  resultant  precision 
fire  artillery  mission. 

Some  of  the  most  important  background  conditions  which 
were  in  effect  when  the  model  was  run  to  obtain  the  vulnera- 
bility data  contained  in  Figures  9,  10  and  11  are  listed 
here . 

AO  dimension  = 3000  meters 

by  3000  meters 

Nap-of -the-earth  maximum  = 45° 

angular  deviation  from 
beeline  ground  track 
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T1  = 20  seconds  appears  essentially  as  a level  line  is 
probably  a function  of  the  fact  that  the  time  between 
artillery  impacts  is  now  getting  very  large  when  compared 
to  increasing  amounts  of  time  being  spent  at  the  initial 
battle  position.  Consequently,  AH  vulnerability  is  becoming 
less  sensitive  to  changing  levels  of  the  CALL. TO. FIRE 
parameter  at  the  higher  T1  values. 

The  results  of  Figure  10  differ  from  the  results  of 
Figure  9 in  the  increased  vulnerability  at  corresponding 
states  of  the  environment.  This  is  because  Figure  10 
contains  output  from  those  replications  of  the  model  in 
which  a forward  observer  and  precision  fire  missions  were 
simulated  in  addition  to  the  area  fire  missions.  Another 
difference  in  the  results  is  that  the  rate  of  vulnerability 
increase  along  the  entire  range  of  CALL.TO.nIRE  values 
changes  dramatically  with  the  different  values  of  Tl. 

At  state  Tl  = 2 seconds,  the  vulnerability  increase  is 
37.1%;  at  state  Tl  = 11  seconds,  the  increase  is  65.4%; 
at  state  Tl  = 20  seconds,  the  increase  is  70.1%.  This  is 
explained  by  the  fact  that,  as  the  intensity  of  the  area 
fire  diminishes  (Tl  getting  larger),  the  threat  to  the  AH 
by  the  precision  fire  missions  becomes  a greater  part  of 
the  total  threat.  Consequently,  longer  stays  in  the  AO, 
which  is  indicative  of  increasing  CALL. TO. FIRE  values,  means 
that  the  overall  AH  vulnerability  becomes  more  sensitive 
to  the  precision  fire  encounters  than  to  the  area  fire 
encounters.  The  increased  rate  of  vulnerability  shown  by 
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Figure  10.  AH  Vulnerability  to  Hypothetical  Artillery 
Threat,  Area  Fire  Plus  Precision  Fire  with 
C3. DELAY  = 60  Seconds  and  TLE  = 50  Meters. 


the  model  is  proof  of  that. 

The  results  of  Figure  11  are  more  difficult  to  explain. 
The  background  conditions  for  this  figure  included  area 
fire  at  five-second  intervals  while  the  CALL. TO. FIRE  value 
was  held  constant  at  50  seconds.  As  shown  by  the  TLE  = 10 
meters  and  TLE  = 100  meters  curves,  the  model  produces  the 
intuitively  correct  response  that  increasing  forward 
observer  accuracy  (TLE  getting  smaller),  the  AH  vulnerability 
increases.  At  once,  the  TLE  = 55  meters  curve  seems  out  of 
place.  One  answer  that  can  be  offered  to  explain  this 
apparent  anomaly  is  that  the  replication  planning  factors 
applied  to  this  research  only  give  the  experimenter  a 90% 
confidence  level  that  the  model's  output  lies  within  ± 0.03 
of  the  actual  for  the  environment  being  simulated.  Since 
the  TLE  = 10  meter  and  TLE  = 55  meter  curves  are  everywhere 
within  a of  0.01  of  each  other,  there  remains  some 
possibility  that,  for  all  practical  purposes,  the  vulner- 
ability results  of  Figure  11  are  insensitive  to  changing 
levels  of  TLE  in  the  range  of  10  meters  to  55  meters. 
Additional  experimentation  with  the  model  is  warranted  in 
order  to  more  completely  understand  the  interaction  among 
variables.  Another  objective  for  further  experiments  might 
be  to  determine  why  two  of  the  curves  in  Figure  11  are  shaped 
concave  upward  while  the  lower  curve  is  concave  downward. 

At  any  rate,  the  results  of  Figure  11  are  intuitively 
correct  regarding  the  C3. DELAY  parameter.  What  the  model 
predicts  is  that  increasing  the  time  delay  between  target 
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Figure  11.  AH  Vulnerability  to  Hypothetical  Artillery 
Threat,  Area  Fire  Plus  Precision  Fire  with 
T1  = 5 Seconds  and  CALL. TO. FIRE  = 50  Seconds 
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detection  and  initiation  of  the  precision  fire  mission 
generally  decreases  the  AH  vulnerability  to  the  artillery 
threat . 

This  chapter  was  certainly  not  intended  to  be  an 
exhaustive  example  of  this  model's  capabilities.  It  was, 
however,  meant  to  show  the  general  form  of  the  model's 
output.  The  user  of  the  model,  by  careful  selection  of  the 
input  parameter  values,  will  be  able  to  examine  in  detail 
as  much  of  the  artillery  versus  AH  environment  that  he 
so  desires. 
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properly.  This  goal  was  met  by  rigorous  examination  of  the 
model's  output  in  graphical  form  and  from  the  exchange  of 
ideas  between  the  researcher  and  experts  in  the  field. 

These  procedures  are  described  in  Chapter  III. 

The  third  goal  was  to  perform  sensitivity  analysis 
upon  the  main  parameters  of  interest.  This  goal  was  met  by 
selecting  four  parameters  and,  then,  experimenting  with  the 
model  while  using  various  levels  of  these  parameters.  The 
model's  results  from  this  experimentation  are  included 
in  Chapter  V. 

The  final  goal  was  to  draw  general  conclusions  about 
the  AH  vulnerability  to  the  artillery  threat  based  upon  the 
model's  output.  This  goal  was  met  by  the  interpretation  of 
the  model's  vulnerability  predictions.  As  was  discussed  in 
Chapter  V,  the  AH  becomes  more  vulnerable  the  longer  it 
remains  operating  in  the  AO.  That  vulnerability  also 
increases  with  increasing  intensity  of  the  simulated  area 
fire.  The  AH  also  becomes  more  vulnerable  as  the  forward 
observer  is  credited  with  having  better  C3  capability 
with  his  supporting  artillery  unit.  The  model  generally 
shows  that  increasing  AH  vulnerability  is  also  a function 
of  increased  forward  observer  accuracy  in  specifying  the 
AH's  location.  All  of  these  predictions  were  as  expected. 
The  fact  that  the  user  of  this  model  has  numerous  variables 
to  specify  prior  to  experimentation  makes  the  forming  of 
definitive  conclusions  about  the  model's  performance 
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difficult  to  make  other  than  in  these  most  general  of 
terms . 

Recommendations 

As  this  research  effort  drew  to  a close,  reflection 
upon  what  was  accomplished  logically  led  to  several  ideas 
which  could  enhance  the  realism  of  this  simulation  model. 
These  ideas  are  described  in  this  section. 

Currently,  the  model  simulates  the  TOW  engagement 
sequence  as  simply  a fixed  amount  of  time  (35  seconds) 
that  the  AH  is  required  to  remain  at  a battle  position 
prior  to  moving  on.  This  treatment  obviously  does  not 
capture  very  much  of  what  actually  transpires  during  this 
AH  activity.  One  idea  for  improving  this  part  of  the  mode 
is  to  allow  the  AH  to  be  characterized  by  some  specified 
height  above  the  ground  during  the  TOW  engagement.  This 
additional  attribute  for  the  AH  is  in  contrast  to  the 
model's  current  handling  of  AH  position  in  two  dimensions 
only.  This  height  parameter  would  then  influence  the 
probability  that  line  of  sight  exists  and  the  probability 
of  detection  by  the  forward  observer.  The  fact  that  a 
TOW  launch  has  or  has  not  occurred  could  become  another 
random  variable  and  the  missile  launch  itself  should  have 
some  bearing  upon  the  AH's  detectability.  Another  idea 
is  to  subject  the  AH  crew  to  similar  line  of  sight  and 
detection  algorithms  in  the  crew's  search  for  the  armor 
target.  The  effects  of  modeling  this  battlefield  activity 
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surely  would  cascade  into  the  TOW  launch  versus  no-launch 
concept  and  into  the  possibility  of  whether  or  not  subsequent 
reattacks  are  to  be  made.  Also,  the  present  model  has 
constrained  the  number  of  TOW  engagements  per  sortie  in  the 
AO  to  a maximum  of  four.  However,  since  the  aircraft  is 
capable  of  being  loaded  with  eight  missiles,  the  maximum 
number  of  TOW  engagements  available  to  the  AH  during  any 
single  replication  should  be  increased  to  eight  with, 
perhaps,  different  probabilities  associated  with  each  weapon 
load . 

At  this  state  of  model  development,  the  movement  of  the 
AH  within  the  AO  is  not  influenced  in  any  way  by  the  artillery 
activity.  A logical  improvement  to  this  is  to  incorporate 
decision  rules  for  the  simulated  AH  crew  to  follow  as  the 
crew  is  given  the  capability  to  sample  the  artillery  environ- 
ment and  react  to  it.  These  decision  rules  could  be  flexible 
to  the  extent  that  the  decision  logic  is  adaptive  to  the 
urgency  of  the  AH's  anti-armor  mission.  Candidate  cues  for 
the  AH  crew  to  use  in  exercising  their  newly  created  decision 
making  capability  are  the  proximity  of  detonating  rounds,  the 
frequency  with  which  rounds  are  detonating  in  the  crew's 
available  field  of  view,  the  aspect  of  the  encounter  (front, 
right  rear  quarter,  etc.),  and  interactions  among  these  three 
cues.  The  inclusion  of  aspect  would  necessarily  require  an 
upgrade  to  the  AH's  positional  descriptor  from  two  dimension 
coordinate,  as  is  presently  simulated,  to  two  dimension 
coordinate  plus  an  azimuth  reference  for  the  aircraft's 
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longitudinal  axis.  The  goals  of  such  a modification  should 
be  to  prevent  the  AH  entry  into  a region  which  is  experiencing 
a heavy  artillery  preparation  and  to  command  a hasty  departure 
from  such  a region  should  the  AH  be  caught  unawares. 

At  present,  only  a single  gun  is  modeled  during  the 
precision  fire  artillery  missions.  A suggestion  offered  by 
the  Fighter  Division,  Directorate  of  Studies  and  Analyses, 
Headquarters  Air  Force  is  to  upgrade  this  one  gun  attack  to  a 
multiple  gun  attack  by  using  the  concept  of  Mean  Point  of 
Impact  (MPI).  This  methodology  essentially  takes  into  account 
that  the  central  aimpoint  of  a single  weapon,  from  a battery 
of  weapons,  will  be  randomly  distributed  about  an  MPI  in  a 
manner  analogous  to  that  which  is  presently  used  in  the 
Event  THREAT  by  probable  errors  in  range  and  deflection.  After 
the  central  aimpoint  for  each  weapon  has  been  specified,  the 
application  of  individual  rounds  should  proceed  as  presently 
formulated . 

One  final  recommendation  is  for  an  examination  of  other 
measures  of  AH  vulnerability  that  may  be  just  as  appropriate 
for  the  model's  user  as  that  which  is  contained  in  this  model. 
This  model  specifies  AH  vulnerability  to  be  the  computed  P^ 
per  AH  sortie  within  the  AO.  Other  choices  for  vulnerability 
may  be  P^  per  unit  time  and  P^  per  number  of  TOW  engagements 
which  were  attempted.  Of  the  four  general  recommendations 
presented  in  this  section,  this  last  recommendation  appears 
to  weigh  most  heavily  upon  the  particular  user's  preference. 
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APPENDIX  A 


HYPOTHETICAL  THREAT  DATA  AND  CODE  LISTING 

The  first  page  of  this  appendix  contains  the  vulner- 
ability curve  of  the  hypothetical  artillery  weapon  as  a 
function  of  range  between  AH  and  burst  point  of  a round. 
The  remaining  pages  of  this  appendix  contain  the  code 
listing  for  this  computer  simulation  model. 
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JOB  CARO 
CONTROL  CAROS 
ENO  OF  RECORD  CARO 


SUBROUTINE  SCRinE(AHX,AHY,ARTYX,ARTYY,XMAX,YMAX,K,L,J,U,PKX,PKY, 
CPKS,NRUNS,PER,PFD,AIMX,AIMY,IFLAG> 

DIMENSION  AHX(l) , AHY<1) , ARTYX(l), ARTYYI1) ,PKX<1) ,PKY(1) ,PKS<1> 

XC  *XMAX/4. 5 
YC=YMAX/4.5 

AHX(K*1)=0.  $AHX<K*2)=XC 
AHY (K*i) =0 • JAHY(K*2)=YC 
ARTYX  <L*i) =0 • $ARTYXtL*2)*XC 
ARTYY  (L  *1)  =0  , f ARTYY  (l  *2)  = YC 
PKX<U+1)  = 0.  $PKX<I1*2> = XC 

PKY(I1+1>=0.  IPKY  <I1*2)  =YC 
F=  FLOAT  <JI 

CALL  PLOT  ( 0 • , -3  . , -3) 

CALL  PLOT<2.,2.5,-3> 

60  CALL  SPAXIS(0.,0.,15HX-AXIS  ( METE  RS  > , -IF  , «» . 5 , 0 . , 0 . , XC , 1 . IT  , -.5,  . 15 
C|0«t iTSi Q)  0«) 

CALL  SPAXIS<0. ,D.,15HY-AXIS  ( METERS  ) , 15  ,4  .5 , 90.  , 0.  , YC,-  .687 , 1 .17, 
C. 15, 90. ,.75, 0,0.) 

CALL  PLOT<0.,4.5,3> 

CALL  PL0T<4.5,4.5,2) 

CALL  PLOT (4. 5*0. *2) 

CALL  PL0T(4. 9625,-1. ,3) 

CALL  PLOT(-1.12c,-l.,2) 

CALL  PLOT( -1.125,5,5,2) 

CALL  PL0T(4.962r,5.5,2) 

CALL  PL0T<4. 9625,-1. ,2) 

CALL  PLOT<4.9725,-1.01,3) 

CALL  PLOT(-1.136,-1.0i»2) 

CALL  PLOT(-l. 136,5.51, 2) 

CALL  PLOT (4. 9725, 5. 51, 2) 

CALL  PLOT(4.9725,-1.01,2) 

CALL  SYMBOL ( 1.1 25, 4. 625,. 375, 3HRUN, 0.,3> 

CALL  SYMBOL (2.5,4.625, .375,75, 0.,-l> 

CALL  NUMBER (3. 0,4. 625,. 375, F,0. ,-ll 
IF<L40.EQ.8>70,R0 
60  CALL  LINE(  AHX,  AHY,  K,  1,  0, 11) 

IF  III.GE.I) 10, 2P 
10  DO  100  1=1,11 

XI =PKX 11)7 XC 4.15 
Y1=PXY  (D/YC-.07 
X2=PKX < I ) /XC 
Y2*PKY ( I ) FYC 

CALL  SYHBOHX2,Y2,.07,0,0.,-1) 

CALL  SYMBOL(X1,Y1,.07,2HPX,0.,2) 

CALL  SYMBOL  <999.,  999. , . 07 ,7 1, 0 . ,-l) 

100  CALL  NUM9ER<999.,999. , .07.PKS <I),0. ,6) 

20  CALL  PLOT (11., 0., -3) 

L4  0=8 
GO  TO  60 

70  CALL  LINE(ARTYX,ARTYY,L,1,-1,3) 

IFdFLAG.EO. 1)71, 40 

71  Al=(AIMX-1.5B9*r£0)/XC 
A2  = <AIMY  *1 .559* PER  I /YC 
B1  = <AINXM.559*PE0)/XC 
82*A2 

C1«B1 

C2* < AI NY -1.559* PER )/YC 
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01*Ai 

D2=C2 

CALL  PL0T(A1,A2,3> 
CALL  PL0T(B1,B2,2) 
CALL  PLOT(Ci,C2,2> 
CALL  PLOT (01,02,2) 
CALL  PLOT  ( Aj  ,A2,2) 

AO  CALL  PLOT(il.,0.,-3) 
95  ENO 


END  OP  RECORO  CARO 


PREAMBLE 

LAST  COLUMN  IS  BO 

DEFINE  I , J, K,L, A, II, NO. RUNS, 12, FLAG, FLA G1 ,ELAG2, FT, OPTION 
AS  INTEGER  VARIABLES 

OEFINE  X.MAX,Y.MAX,BP.X,8P.Y,ENTRYy , ENT RY V, XSLOHOOWN, YSLOHDOHN, XSPEEOUP, COUNT, 
YSPEEDUP.EXIT.X, EXIT.Y, 1AITX , 1ALTY,  2 ALT  X, 2ALTY, 3ALTX, SALT Y, PK1, AIMX, AIMY, 

HI SPEED, LOSPEEO, OEL .T , PHI. MAX, CALL. TO. FITE, XI , Yi , D. X, D. Y, THET A. I , R. I ,T1 , PK. STEP, 
C3. OEL AY, PE.R, PE. D, TL E, ROF,ROUNOS,FOX ,FOf , TERRA IN,H, X22 , Y22 , RNG , AL 1,01, 

FO.YLINE  AS  REAL  VARIABLES 

DEFINE  AHX, AMY, ARTYX.ARTYY.PK. THREAT, PK.X,PK.Y,PKS,RC  AS  i-DIM  ARRAYS 

DEFINE  IL  AS  AN  INTEGER, 1-OIM  ARRAY 

OEFINE  SCRIBE  AS  A FORTRAN  ROUTINE 

EVENT  NOTICES  INCLUDE  NOE, SCENARIO, I NCR  ESS, THPEAT .ENTRY . L EG , 1 .RE ATTACK, 

2. RE ATTACK, 3. REATT ACK, DROPBACK, EGRESS, SI  MEND, PREC. FIRE, LOS, DETECT 
TEMPORARY  ENTITIES 

EVERY  AH  HAS  SOME  X.SOME  Y,  A TIME,  AND  BELONGS  TO  A POSITION 
THE  SYSTEM  OHNS  THE  POSITION 

TALLY  ST  AT. 1 AS  THE  HFAN  AND  STAT.2  AS  THE  VARIANCE  OF  PK1 
END 


MAIN 

RESERVE  ARTYX  (*)  ,AND  ARTYY  (*)  AS  999, AND  PK,X  ( * ) , PK.  Y (•  ) , AND  PKS(’)  AS  999 
RESERVE  IL  (* ) , AND  RC(»)  As'  2 

••INPUT  USER-DESIRED  VARIABLES 

READ  X. H AX, Y. MAX, DEL. T, PHI. MAX, Tl, A, PK. STEP, HISPEEO, LOSPEEO, CALL. TO. FIRE, NO. RUNS 
READ  C3.DELAY,PE.R,PF.D,TLE, ROF , ROUNDS, FO . YLINE , T ERRAIN ,H 
REAO  AL1, OPTION 

••X. MAX, Y. MAX  ARE  THE  AO  OTMENSIONS  IN  METERS 

•*OEL.T  IS  THE  TIME  INCREMENT  FOR  NOE  VELOCITY  VECTOR  IN  SECOMOS 
••PHl.MAX  IS  THE  MAX  ANGULAR  OEFLECTION  OF  THE  VELOCITY  VECTOR  IN  DEGREES 
••Tl  IS  THE  ARTY  FIRING  INTERVAL  IN  SFCONOS 
••A  IS  THE  NUM9ER  OF  PK  DATA  POINTS  TO  BE  ENTERED 

••PK.STEP  IS  THE  PANGE  BETWEEN  EACH  SUCCESSIVE  PK  DATA  POINT  IN  METERS 
••NO. RUNS  IS  THE  DESIRED  NUMPER  OF  REPLICATIONS  FOR  THIS  SINGLE  JOB 
••HISPEEO  IS  THE  VELOCITY  USED  DURING  I NGR ESS/ EGRESS  IN  KNOTS 
••LOSPEEO  IS  THE  UPPER  BOUND  VELOCI IY  FOR  NOE  FLIGHT  IN  KNOTS 
••CALL. TO. FIRE  IS  THE  TIME  DELAY  AFTER  ARRIVAL  AT  THE  INITIAL  BATTLE 
••POSITION  PPIOR  TO  THE  UNMASK  MANEUVER  IN  SECONDS 
••C3.0ELAY  IS  TIME  LOST,  DETECTION  T3  1ST  ROUND,  IN  HINUTES 
••PE.R  IS  PFOBABLF  ERROR  IN  RANGE  (METERS) 

••PE.O  IS  PPOBABLE  ERROR  IN  OFELECTIDN  (METEP.S) 

•*TLE  IS  TARGET  LOCATION  ERROR  IN  M FT  EPS 

••ROE  IS  SINGLE  GUN  RATE  0*"  FIRE  IN  ROUNOS  PER  HINUTE 

♦•ROUNDS  IS  TOTAL  AMOUNT  OF  POUNOS  TJ  BE  EXPENDED  OURING  PRECISION  EIRE 
••FO.YLINE  IS  DISTANCE  ABOVF  TOP  OF  AO  OF  THE  FO  POSITION  IN  METERS 
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••TERRAIN  I "5  THE  TERRAIN  TYPE  COEFFICIENT 

••H  IS  THE  FO'S  VIEWING  HETGHT  AGL  IN  MET FRS 

••All  IS  THr  TARGFT  ILLUMINATION  IN  FOOT-CANOLFS 

••IF  OPTIONS  PREC  FIRE  IS  CONDUCTED.  OPTION=0  APEA  FIRE  ONLY. 

PRINT  9 LINES  WITH  X. HAX ,V .MAX, DEL . T , PHI . MA X, T1 , HISPEED,LOSPEEO, CALL . TO . FIRE, 
NO. RUNS  AS  FOLLOWS 

#*«*#  USER  HAS  INPUT  THE  FOLLOWING  PARAMETERS  MAM 
AO  DIMENSIONS  ARE  *»*».  METERS  BY  *•*».  METERS 
NOE  DISPLACEMENT  STEPS  ARE  TAKEN  AT  **. -SECOND  INTERVALS 
MAX  ANGULAR  DEVIATION  FROM  BEE-LINE  PATH  IS  •*.  OEGREES 
ARTILLERY  IS  IMPACTING  AT  ".-SECOND  INTERVALS 
INGRESS/EGRESS  AH  VELOCITY  IS  KNOTS 

NOE  HAX  SPEED  LIMITEO  TO  *».  KNOTS 

AT  INITIAL  PP,  •**.  SECONDS  OF  DELAY  IS  DIRECTED  PRIOR  TO  UNMASK 
**».  REPLICATIONS  ARE  OESIREO 

PPINT  9 LINES  WITH  C3.0ELA Y , PE .R, PE . 0 , TLE ,ROF .ROUNDS , FO.YLINE , TERRAIN, H AS 
FOLLOWS 

C-CU9EO  DELAY  IS  *.*  MINUTES 

PROBABLE  ERROR  IN  RANGE  IS  **•.  METERS 

PROBABLE  ERROR  IN  DEFLECTION  IS  ••*.  METERS 

TARGET  LOCATION  ERROR  IS  METERS 

SINGLE  WEAPON  RATE  OF  FIRE  IS  ROUNDS  PER  MINUTE 

TOTAL  ROUNDS  PROVIOEO  FOR  THIS  MISSION  IS  *". 

FO  IS  RANDOMLY  LOCATED  ALONG  A LINE  FROM  THE  TOP  OF  THE  AO 

THE  TERRAIN  TYPE  COEFFICIENT  IS  *,* 

THE  OBSERVER'S  ALTITUDE  IS  •• . METERS  AGL 
PRINT  2 LINES  WITH  AL1* OPTION  AS  FOLLOWS 

THE  ILLUMINATION  AT  THE  TGT  IS  »***.*»**♦  FOOT-CANDLES 
THE  PRECISION  FIRE  OPTION  IS  • 

RESERVE  RC.IL  AS  2 
RESERVE  PK. THREAT  AS  A 
READ  PK. THREAT 
REAO  FT 
IF  FT  EO  1, 

LET  PK.STEP=PK.STEP* 0.3048  "CONVERTS  FNTEREO  DATA  FROH  FEET  TO  METERS. 
ALWAYS 

LET  DEL.T=DEL.mo. 

LET  T 1=T 1/60. 

LET  LOSPEED=LOSPEEO* .C14A4 
LET  HISP EEO=HI SPEED*  • c 14<*4 
LET  PHI . MAX=PHI .HAX/RAOIAN.C 
LET  J=i 

SCHEDULE  A SCENARIO  NOW 

SCHEOULE  A THREAT  IN  UNIFORM. F ( 0 . 0, T 1 ,5 ) MINUTES 
IF  OPTION  EO  1, 

SCHEOULE  AN  LOS  IN  1./6.  MINUTES 
ALWAYS 

START  SIMULATION 
ENO 


EVENT  SCENARIO 
LET  1=1*1 

IF  1=1, SCHEOULE  AN  INGRESS  NOW 
ALWAYS 

IF  1=2, SCHEOULE  A 1 .RCATTACK  IN  (CALL. TO. FI  PE *35. )/60.  MINUTES 
ALWAYS 

IF  1=3, SCHEOULE  A ?. REATTACK  IN  35. /60.  NINUTFS 
ALWAYS 

IF  !=«,, SCHEDULE  A 3. REATTACK  IN  35.^60.  MINUTES 
ALWAYS 

IF  1 = 5 ANO  12=0, SCHEDULE  A ORQPRACK  IN  <5  ALL . TO . FIRE ♦ 35 . > /60 . MINUTES 
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always 

IF  1 = 5 AND  12  HE  0,SCHEOULE  A OROPBACK  IN  35. /60.  MINUTES 
ALWAYS 

IF  I=6tSCME0ULE  AN  EGRESS  IN  R.:FHISPEEOF60.  MINUTES 

ALWAYS 

RETURN 

END 


EVENT  INGRESS  "TO  DETERMINE  THE  PARTICULAR  SCENARIO  OF  THIS  RUN 
LET  BP.X=UNIFORM.F(D.0,X.MAX,l> 

LET  BP.Y=UNIFORM.F(Y.MAX-m  .,Y.HAX,i> 

LET  ENTRYX=UNI FORM. FCO.O.X. MAX, 1) 

LET  ENTRYY=0. D 

LET  FOX= UNI  FORM.  F (0.0,X.MAX,1) 

LET  FOY=Y.HAX*FO.YLINE 
LET  YSLOWDOWN=Y. MAX-1000. 

••FROM  SIMILAR  TRIANGLES.... 

LET  XSLO WOOWN=ENTR YX- ( ENTRY X -BP. X) • YSLO WOOWNFBP . Y 

••THROW  THE  OICE  TO  DETERMINE  THE  NUMBER  OF  REATTACKS  TO  ACCOMPLISH,  IF  ANY. 
LET  YSPEEOUP= YSL OWDOWN 
LET  EXIT  • > =0  • 0 
LET  I 2=RANDI.F (0,3,2) 

IF  12=0, LET  1=4 
LET  XSPEEDUP=BP.X 
LET  EXIT  .X=BP , X 
ALWAYS 

••BIAS  ALL  SUBSEOUENT  RE-ATTACKS,  IF  ANY,  TO  THE  CENTER  OF  THE  AO. 

IF  BP.X  LT  X.MAX/2,,LET  BIAS=1. 

ELSE 

LET  BIAS*-1. 

ALWAYS 

••SUBSEQUENT  RE-ATTACKS  ARE  ALLOWED  TO  VARY  FROM  THE  PREVIOUS  ATTACKING 
••POSITION  BY  ♦ OR  - 50  METERS  IN  THE  Y COORDINATE  AND  UNIFORHLY 
••VARYING  FPOM  30  METERS  TO  150  METERS  IN  LATERAL  DISPLACEMENT. 

IF  12=1, LET  1ALTY=UNIFORM.FIBP.Y-5J).,BP.Y+50.,1) 

IF  1ALTY  GT  Y. MAX, LET  1ALT Y=Y.HAX 
ALWAYS 

LET  1 ALTX=BP. X tBI AS* UNIFORM. F( 30., 150.,  1) 

LET  XSPEE0UP=1ALTX 
LET  EXIT  »X=1ALTX 
ALWAYS 

IF  12=2, LET  1ALTY=UNIFORM.F(BP.Y-50.,BP.Y*50.,1) 

IF  1ALTY  GT  Y.MAX.LFT  1ALTY=Y.MAX 
ALWAYS 

LET  1ALTX=BP.X*BIAS*UNIFORM.F(30.,150.,1» 

LET  2ALTY=UNIFORM.F  tlALTY-50 . , 1ALTY *50. , 1 ) 

IF  2ALTY  GT  Y. MAX, LET  2ALTY=Y.MAX 
ALWAYS 

LET  2ALTX=1ALTX«PIAS*UNIFORM.FC3P.,150.,1> 

LET  XSPEEOUP=2ALTX 
LET  EXIT .X=2ALTX 
ALWAYS 

IF  12=3, LET  1ALTY=UNIFORM.F(BP.Y-50.,BP.Y*50.,1) 

IF  1ALTY  GT  Y. MAX, LET  1ALTV=Y.MAX 
ALWAYS 

LET  1 ALT  X=BP.  X »R  IAS*  UNIFORM.  F(  TO.,  15  0.,  1) 

LET  2ALTY=UNIFOPN.FUALTY-50.,1ALTY*50.  ,1) 

IF  2ALTY  GT  Y. MAX, LET  2ALTY= Y. MAX 
ALWAYS 

LET  2ALT X=1ALTX» PI  AS* UNIFORM. F( 30.,  15  0. ,1) 

LET  3ALT Y*UNIFORM. F (2ALTY-50.,2ALTY»50.,U 
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IF  3ALTY  GT  Y. MAX, LET  3ALTV=Y.MAX 
ALWAYS 

LET  3ALTX*2ALTX*BIAS*UNIFORM.FI30.,15O.,l> 

LET  XSP£E0UP*3Al  TX 
LET  EXIT.X*3ALTX 
ALWAYS 

LET  O.X*XSLOWDOWN 
LET  D.Y*YSLOWDOWN 

•'INITIALIZE  THE  AH  POSITION  AT  THE  ENTRY  POINT  TO  THE  AO 
CREATE  AN  AH 
LET  X*ENTRYX 
LET  Y=0.0 
LET  X1*X 
LET  Yi*Y 

LET  TIME=TIME.V'1440. 

FILE  THIS  AH  IN  POSITION 

LET  R.I*SaRT.FHENTRYX-XSL0W00NN>**2*YSL0WD0WN*»2) 

LET  THETA. 1*0,0 

SCHFOOLE  AN  ENTRY. LFG  IN  R.I /HISPEED/60 . MINUTES 

RETURN 

END 


EVENT  ENTRY. LEG  “EMPLOYS  STRAIGHT  LINE  FLIGHT 

CREATE  AN  AH 

LET  X*D. X 

LET  Y=O.Y 

LET  X1*X 

LET  Y1*Y 

LET  TIHE*TIM£«V*1440« 

FILE  THIS  AH  IN  POSITION 
LET  D.X=BP.X 
LET  D.Y=BP.Y 

LET  R.I=OEL.T*60. •UNIFORM. F ( ,5'LOSPEEO, LOSPEEO, 4) 

LET  THETA.I*UNIFORM.F(-l.»PHI.MAX,PHI.NAX,4> 

SCHEDULE  AN  NOE  NOW 

RETURN 

END 


EVENT  1. REATTACK 
LET  O.X*  1ALTX 
LET  O.Y*  1ALTY 

LET  R. I*OEL. T*60. ’UNIFORM. F ( ,5*LOSPEEO, LOSPEEO, 4) 

IF  12*1, LET  1*4 

ALWAYS 

SCHEDULE  AN  NOE  NOW 

RETURN 

END 


EVENT  2. REATTACK 
LET  O.X* 2ALTX 
LET  O.Y*2ALTY 

LET  R. I*OEL.T*6n. ‘UNIFORM. F<.5‘tOSPEED,LOSPFED, 4) 

IF  12*2, LET  1*4 

ALWAYS 

SCHFOULE  AN  NOE  NOW 

RETURN 

END 


64 


EVENT  3. REATTACK 
IET  0.X-3ALTX 
LET  C.Y=3ALTY 

LET  R. I=OEL.T*SO. •UNIFORM. F( .5*LOSPEEO, LOSPEFO, 4> 

SCHEDULE  AN  NOE  NON 

RETURN 

END 


EVENT  DROPBACK 
LET  O.X=XSPEEOUP 
LET  D • Y=  YSPEEDUP 

LET  R.I*OEL.T*60. 'UNIFORM. F(.5*LOSPEEO,LOSPEEO, 4) 

SCHEDULE  AN  NOE  NON 

RETURN 

ENO 


EVENT  EGRESS  • ‘EMPLOYS  STRAIGHT  LINE  FLIGHT 
CREATE  AN  AH 
LET  X=D.X 
LET  Y*D.Y 

LET  TIME*TIMS.V*1440. 

FILE  THIS  AH  IN  POSITION 
SCHEDULE  A SIMEND  NOH 
RETURN 
ENO 


EVENT  NOE 

••SIMULATES  NAP  OF  THE  EARTH  FLIGHT  OY  GENERATING  A SUCCESSION  OF  SHORT, 
••STRAIGHT  POSITION  VECTORS  OF  VARYING  LENGTH  ANO  VARYING  ANGULAR 
••DISPLACEMENT  FROM  THE  BEE-LINE  ORIENTATION. 

CREATE  AN  AH 

LET  TIME=TIME.V»1440. 

LET  RANGE=SORT.FUO.X-X1)**2m*».Y-Y1)"*2> 

IF  RANGE  LE  R.I,  LET  X«O.X 
LET  Y*O.V 
LET  X1=X 
LET  V1*Y 

FILE  THIS  AH  IN  POSITION 
LET  R. 1=0.0 

IF  1=5, LET  O.X=XSPEEDUP 

LET  O.Y=0.0 

LET  R.I=YSPEEDUP 

LET  THETA. 1*0.0 

ALHAYS 

SCHEDULE  A SCENARIO  NOH 

GO  TO  CHECK 

ELSE 

LET  X1«X1MR.I*C0S.F<THETA.I)*(D.X-X1»»R.I*SIN.F(THETA.I>*(V1-D.Y>) /RANGE 
LET  Y1«71MR.I*C0S.F (THETA .1)*  (O.Y-Yl)'R. I'SIN.F (THETA .1) * (0. X-Xi) ) /RANGE 
•'THESE  INSTRUCTIONS  INSURE  THAT  THE  AH  REMAINS  HITHIN  THE  AO. 

IF  XI  LT  0.0, LET  X1=0.0 
ALHAYS 

IF  XI  GT  X. MAX, LET  Xi=X.MAX 
ALHAYS 

IF  Y 1 LT  0.0, LET  Y1*0.0 
ALHAYS 

IF  VI  GT  Y. MAX, LET  V1=V.HAX 

ALHAYS 

LET  X*X1 
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LET  Y*Y1 

FILF  THIS  AH  IN  POSITION 

LET  R.I*0EL.T»6P.*UNIF0RH.F(  .5*L0SPEE0,LDSPEE0,  A) 
LET  THETA.I=»UNIFORH.F  <-l  ,»PHI .MAX, PHI .MAX , A) 
SCHEDULE  AH  HOE  IN  OEL.T  MINUTES 
•CHECK*  RETURN 
END 


EVENT  LOS  '‘THIS  EVENT  IS  SCHEOULEO  AT  10  SECOND  INTERVALS  THROUGHOUT  THE 
••LENGTH  or  THE  SIMULATION.  THE  ACTIVITIES  HITHIN  THIS  EVENT 
••DETERMINE  WHETHER  OR  NOT  LI NE-OF-SI GHT  EXISTS. 

IF  R.I  EO  0.0, LET  X2  = X1 
LET  Y2=Y1 
GO  TO  CHECK 
ELSF 

IF  R.I  GT  OEL.T'60.*LOSPEFD,LET  D=R, I/HISPEED/60 . 

ELSE 

LET  0=DEL«T 
ALWAYS 

••THIS  NEXT  SECTION  PINPOINTS  THE  AH  LOCATION  WHEN  THE  AH  IS  IN  HOTION  ANO 
••ARTY  ROUNDS  IHPACT  DURING  THE  EXECJTION  OF  THE  AH  DISCRETE  DISPLACEMENT 
••STEP. 

FOR  EACH  AH  IN  POSITION  IN  REVERSE  ORDER,  FIND  THE  FIRST  CASE 

LET  PART .HOVE= (TIME. V' 1440 .-T I HE) /D 

LET  RANGE=>SORT.FUO.X-Xi>**2«-rD.Y-Yl>*»2) 

IF  RANGE  LE  R.I,  LET  X2*Xl*CO.X-Xi)  *P ART.  HOVE 
LET  Y2=Y1*(0.Y-Y1) »PAPT.MOVE 
GO  TO  CHECK 
ELSE 

LET  X3«X1MR.I*C0S.F<THETA.IJ*  (0.  X-XH  »R.  I*  SIN.F1THET  A.H  *1  Yl-O.  YM  /RANGE 
LET  Y 3*Y  1MR.I  *COS  »F  (THETA, II*  CO.Y-Y1HP.  I*SIN.  F ITHET  A.  I>*  TO.  X-Xl)) /RANGE 
IF  X3  LT  0.0, LET  X3=0.0 
ALWAYS 

IF  X3  GT  X. MAX, LET  X3=X.HAX 
ALWAYS 

IF  V3  LT  0.0, LET  Y3«0.0 
ALWAYS 

IF  Y3  GT  Y. HAX, LET  T3*Y.MAX 
ALWAYS 

LET  X2»X1*U3-XU»OART.MOVE 
LET  Y2=Y1MY3-Y1)*PART.M0VE 

•CHECK*  LET  RANGE».001*SQRT.F< «FOX-X2>**2  MFOY-V2)**  2) 

LET  RBARsTERRAIN-  (1 .56*3  .A6»  C TH/l  00  . ) *•  1.  17  3)) 

LET  PLOS*  C2.*RANGE/RBAR  + 1.)*  EXP.F  C*2. ’RANGE/RBAR) 

LET  OICE=UNIFORM.FTO. 0,1. 0,10) 

IF  OICE  GT  PL OS, SCHEOULE  AN  LOS  IN  1./6.  MINUTES 
ELSE 

LET  X22*X2 

LET  Y22=V2 

LET  RNG=RANGE 

SCHEOULE  A DETECT  NON 

ALWAYS 

RETURN 

(NO 


EVENT  OETECT 
OEFINE  II, LIGHT 
LET  AJ0B*2.5 
LET  ATTN*0.558 
LET  DIM=A. 12 


••ONLY  ENTERED  IF  LINE-OF -SIGHT  EXISTS. 

AS  INTEGER  VARIABLES 

••THIS  ATMOSPHERIC  ATTENUATION  COEFF  IS  FOR  7 KM  VISIBILITY. 


66 


I 


LET  ACON«0.3  ••THIS  IS  THE  TGT /BACKGROUND  CONTRAST  RATIO. 

••THE  ACON  VALUE • TO  A GREAT  EXTENT,  DRIVES  THE  RESULTS  OF  THIS 
• 'ALGORITHM.  INCREASING  ACON  INCREASES  PROB  OF  OETECTION. 

LET  TBARslOOOOO. 

LET  AL=AL1 

LET  SOG=3.0  * *THIS  PARAMETER  IS  THE  SKY  BRIGHTNESS/GROUND  BRIGHTNESS  PATIO 

LET  R=RNG*100D. 

LET  C*ACON/(l.*SOG*  (EXP  ,F(ATTN*R/10  0 0 .1  -l.ll 

LET  AL*AL*0.7 

IF  AL  LT  100., GO  TO  A30 

ELSE 

LET  IL1*1 
LET  IL2=1 
LET  ACK= 100 . 

LET  AL-IOO. 

GO  TO  ATO 

•A30 • IF  AL  GT  .0001, GO  TO  AAfl 
ELSE 

LET  PAJOB=0.0 
LET  TBAR.100000. 

GO  TO  A600 

•A40 • LET  ACKslOO. 

FOP  11*1  TO  6,00 
LET  IL1=I1 

IF  AL  LE  ACK  AND  AL  GF  ACK/10.,GO  TO  A60 
ELSE 

LET  ACK*ACK/10. 

LOOP 

•A60*  LET  IL2=IL1*1 
•A70*  LET  IL  ( 11  =IL  1 

LET  IL(2)=IL2  - •> 

LET  RC (1) =0 .0 
LET  RC(Z>*0.0 
FOR  11*1  TO  Z,00 
LET  LIGHT*IL  CI1) 

IF  LIGHT  EO  1 , GO  TO  BIO 
ELSE 

IF  LIGHT  EQ  Z,GO  TO  B40 
ELSE 

IF  LIGHT  EO  3, GO  TO  B70 
ELSE 

IF  LIGHT  EO  4, GO  TO  B90 
ELSE 

IF  LIGHT  EO  5, GO  TO  B110 
ELSE 

IF  LIGHT  EQ  6, GO  TO  B130 
ELSE 

IF  LIGHT  EO  7, GO  TO  B150 
ELSE 

•BIO*  LET  RC  (ID *Z.7i 
IF  C GT  0.8, GO  TO  B170 
ELSE 

IF  C IE  0.25, GO  TO  B20 
ELSF 

LET  RC (I 1»*2. 8039221* (C**0 .2291015) 

GO  TO  B170 

•B20 • IF  C LT  0.025, GO  TO  B30 
ELSE 

LET  RC(Il)«2.33f 5166-0.0509533/C 
GO  TO  B170 

•B3C  * LET  RCCI1>*0.0 
GO  TO  B170 
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*040  * LET  RC(I1>*2.74 
IF  C GT  0.8, GO  TO  B170 
ELSE 

IF  C LT  0.35, GO  TO  850 
ELSE 

LET  RC(I1W.8667792*(C**0. 2251284) 

GO  TO  0170 

•050*  IF  C LT  0.025, GO  TO  060 
ELSE 

LET  RC(I1)=2.3 26 29 99-8.0 51 4585 /C 
GO  TO  0170 

•060*  LET  RC(I1)=0.0 
GO  TO  8170 

•070*  LET  RC<I1)=2.29 
IF  C GT  0.7, GO  TO  8170 
ELSE 

IF  C LT  0.03, GO  TO  880 
ELSE 

LET  RCTI 1) *C/ ( 0. 04976 05*0. 3579996*0 
GO  TO  8170 

•080*  LET  RC(IH*0.0 
GO  TO  0170 

•090*  LET  RC(I1)=1.8 
IF  C GT  0.7, GO  TO  8170 
ELSF 

IF  C LT  0.05, GO  TO  B100 
ELSE 

LET  RC (Il)=1.6 1493 45-0. 06510 33 /C 
GO  TO  0170 

•0100*  LET  RC (111=0.0 
GO  TO  0170 

•B110  * LET  RC (I 1) =1 . 32 
IF  C GT  0.8,GO  TO  8170 
ELSE 

IF  C LT  0.084, GO  TO  0120 
ELSE 

LET  RC(I1>«1. 336422-0. 1123916/C 
GO  TO  0170 

•0120*  LFT  RC(I 1) =0.0 
GO  TO  0170 

•0130'  LET  RC (1 1)  =0 .<>7 
IF  C GT  0.6,GO  TO  0170 
ELSE 

IF  C LT  0.18, GO  TO  B1'‘0 
ELSE 

LET  RC(I1>*0. 5686696-0. 0968436/C 
GO  TO  0170 

•0140*  LET  RC(I1)=0.0 
GO  TO  0170 

*0150*  LET  RC(I1)=0.14 
IF  C GT  0.8, GO  TO  0170 
ELSE 

IF  C LE  0.5, GO  TO  8160 
ELSE 

LET  RC(1 1) *0. 29-0. 12/f 
GO  TO  0170 

•8160*  LET  RC(I1»=0.0 
•0170*  LOOP 
LET  RC1*RC(1) 

LET  RC2*RC(?I 

LET  RX=RC2MRC1-PC2>* ( ( AL-ACK710 . I S ( «CK-t CK/1 0. ) ) 
LET  RX»RX*7.0 
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LET  S=D1M/IR71000.) 

LET  RX=RX*S 

IF  RX  LT  0.1, LET  RX=0.0 
LET  PA JO 8=0 .0 
LET  TRAR=100000. 

GO  TO  A600 
ELSE 

LET  XX=(RX-AJ09) /(  0.625*AJ0B) 

IT  XX  GT  -45.0, GO  TO  CIO 
ELSE 

LET  PAJOB=0 • 0 
GO  TO  A500 

•CIO*  IF  XX  LT  45.0, GO  TO  C20 
ELSE 

LET  PAJOB=1.0 

GO  TO  A500 

•C20*  LET  F=1 . 0 

IF  XX  LT  O.O.LET  F=0.0 

ALWAYS 

LET  S2=XXF1.414 
LET  X2=ABS.F<S2) 

LET  E1=1.*<0.0000430638*X2*0.0002765672)»X2»*5 
LET  E1=E1*(0.  000 15 201 34* X2* 0.0 092705 272) *X2»*3 
LET  E1*E1K 0.042 28 201 23*X2* 0.07 052307 84 )*X2 
LET  E1=E1**16 
LET  E1>1.-1./E1 
LET  PAJO8*0.5*El/2. 

IF  F GT  0.0, GO  TO  A500 
ELSE 

LET  PAJOB*l.-PAJOB 

•A500*  LET  TS=1. 7*600. /64. 

LET  PS=1.-EXP.F(-1.7*PX78.) 

IF  PS  GT  0.0, GO  TO  A160 
ELSE 

LET  TBAR=100000. 

GO  TO  A600 

•Alt  0 * LET  TBAR  = 0.5*TS*(2.-PS)/PS 

•A600  * LET  P10=PAJOB*<1.-EXP.F«-10./TBAR)) 

LET  OICE«UNIFORM.F(0.n, 1.0,10) 

IF  MCE  GT  P10, SCHEDULE  AN  LOS  IN  1./6.  MINUTES 
ELSF 

LET  FLAG1*1 
LET  FLAG2*1 

SCHEDULE  A PREC. FIRE  NON 
ALWAYS 

return 

END 


EVENT  PREC. FIRE  "THIS  EVENT  CONTROLS  THE  EMPLOYMENT  OF  A FO-DIRECTEO  FIRE 
••MISSION  AGAINST  A POINT  TARGET. 

IF  FLAG1=1,LET  FLAG1=0 

LET  BETA  *UNIFORM.F  (0.0,2.*PI.C,8) 

LET  PHO=TLE*SORT .F ( UNIFORM.F 10.0,1.0,7) ) 

••FROM  THE  GEOMETRY  OF  THE  PROBLEM 
LET  AIMX«X2?»RHO*COS.F(BETA) 

LET  AIMY=Y22»RH0*SIN.F(BETA) 

LET  FLAG*1 
LET  COUNT* 1 

SCHEDULE  A THREAT  IN  C3.0FLAY  MINUTES 
ELSF 

IF  COUNT  LT  ROUNDS, LET  COUNT=COUNT* 1. 0 
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LET  FLAG*1 

SCHFOULE  A THREAT  TN  l.O/POF  MTNUTES 
ELSE 

LET  COUNT=0.0 

SCHFOULF  AN  LOS  NON 

ALWAYS 

ALWAYS 

RETURN 

END 


EVENT  THREAT 

•'THIS  EVENT  LAYS  OOWN  A UNIFORM  BARRAGE  PATTEPN  WITHIN  THE  AO. 

"IF  REQUIRED,  THIS  EVENT  ALSO  LAYS  DOWN  THE  PRECISION  FIRING  PATTERN 
"IT  ALSO  TFSTS  THE  PPOXIHITY  OF  THE  AH  TO  THE  ARTY  IMPACT  AND  IF  WITHIN 
"THE  IMPACT  THRFAT  ENVELOPF,  CALCULATES  AND  STORES  THE  APPROPRIATE  P-K  DATA 
DEFINE  INDEX  AS  AN  INTEGER  VARIABLE 
LET  L=L ♦ 1 

IF  FLAG  EO  1,LET  FLAG=0 

LET  ART YX (H = AIMX* NORMA L. F ( O.O.PF.D/O.S TV F,F) 

LET  APTYY(L)*AIMY'NORHAL.F(0.0,PF.R/0.6  7I»5,9) 

SCHEOULE  A PREC.FIRE  NOW 

GO  TO  JOB 

ELSF 

LET  ARTYX(L>*UNIFORH.F(0.0,X.HAX,5J 
LET  ARTYY<L)-UNIFORM.F(0.0,Y.MAX,5> 

SCHFOULE  A THREAT  IN  T1  MINUTES 
•JOB* 

"THROUGHOUT,  COORDINATE  (XI, Yl)  ARE  THE  AH  CURRENT  POSITION,  OR  FOR  THE  AH 
"IN  MOTION,  THE  MOST  RECENTLY  OCCUFIEO  POSITION  AHONGST  THE  SERIES  OF 
•'DISCRETE  DISPLACEMENT  STEPS. 

IF  F.I  EO  0.0, LET  X?=X1 
LET  YZ=Y 1 
GO  TO  CHECK 
ELSE 

IF  R.I  GT  OEL.T'60.»LOSPEEO,LET  D=R.I/HlSPEED/60 . 

ELSE 

LET  D*DEL.T 
ALWAYS 

•'THIS  NEXT  SECTION  PINPOINTS  THE  AH  LOCATION  WHEN  THE  AH  IS  IN  MOTION  AND 
••ARTY  ROUNOS  IMPACT  DURING  THE  EXECJTION  OF  THE  AH  DISCRETE  DISPLACEMENT 
•'STEP. 

FOR  EACH  AH  IN  POSITION  IN  REVERSE  ORDER,  FIND  THE  FIRST  CASE 

LET  PART. NOVE=( TIME. V' 1440. -T IMF) /O 

LET  RANGE=SQRT.F((0.X-X1)**2*('I.Y-Y1)"?) 

IF  RANGE  LE  R.I,  LET  X 2*X1» (0. X-Xl) ’PART . MO VE 
LET  Y2«Y1MO.Y-Y1)*PAPT.MOVE 
GO  TO  CHECK 
ELSE 

LET  X3*XiMR.I*COS.F(THETA.I>*tO.X-Xl»'R.I*SIN.F(THETA.I»*(Yl-O.Y) ) /RANGE 
LET  Y3=YlMR.I*C0S.FCfHETA.I>’M0.Y-Yl>*R.I'SIN.F(THETA.I>MD.X-Xin/RANGE 
IF  X3  LT  0.0, LET  X3*0.0 
ALWAYS 

IF  X3  GT  X. MAX, LET  X3=X.MAX 
ALWAYS 

IF  Y3  LT  0.0, LET  Y3*0.0 
ALWAYS 

IF  V3  GT  Y. MAX, LET  Y3=Y.HAX 
ALWAYS 

LET  X2»X1MX3-X1)'PART.M0VE 
LIT  Y?=VI'(Y3-Y1)*®ART.M0VE 

•CHECK'  LET  RANGF  = S0RT.F<(ARTYX<L)-X2»"2MARTYY(L>-Y2)"Z> 
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IF  FANGE  IT  (A-l .) »PK. STEP, LET  INOEX=TRUNC.  F ( RANGE/PK • STEP) 

LET  11=11*1 

LET  C=RANGE-INOEX*PK.STEP 

LET  PK=PK. THREAT < I NCFX*?>  ♦ ( 1 . -C/PK.  ST  EP  > * »PK.T  HRE  AT  (I  NO  EX  *1 > -PK.  T HREA  T ( I NOE  X*?  > ) 
LET  PK.X(I1)=ARTYX(L) 

LET  PK.Y<I1)=ARTYY(H 

LET  PKS(I1)=PK 

ALWAYS 

RETURN 

END 


EVENT  SI  MEND 

OEFINE  G AS  AN  INTEGER  VARIABLE 
IF  II  EO  0, 

LET  PK1=0.0 
ALWAYS 
IF  II  GE  1, 

FOR  G=1  TO  II, 00 

••THIS  DO-LOOP  CHECKS  THE  MONTF  CARLO  TEST  OF  THE  P-K. 

LET  DICEaUNIFORM. Fro. 0,1.0,101 
IF  PKS(G)  GE  OICE, 

LET  PK1=  1.0 
GO  TO  CHECK 
ELSE 

IF  G=I1,LET  PK1=0.0 

ALWAYS 

LOOP 

•CHFCK* 

ALWAYS 

RESERVE  AHX  < • ) , ANO  ANY  <*>  AS  N. POSITION*? 

FOR  EACH  AH  IN  POSITION, DO 

LET  K=K*1 

LET  A HX ( K) =X 

LET  AHY(K»=T 

LOOP 

IF  J«l, 

••WITH  THIS  LOGICAL  IF  STATEMENT,  USER  MAY  SPECIFY  WHICH  RUN  OR  GROUP  OF 
••RUNS  WILL  RE  PLOTTED  WITH  SUBROUTINE  SCRIBE. 

FOR  EACH  AM  IN  POSITION, FINO  THE  FIRST  CASE 
LET  B1=TIME 

FOR  EACH  AH  IN  POSITION  IN  REVERSE  OROEP,  FINO  THE  FIRST  CASE 
LET  B?*TIME 

PRINT  1 LINE  WITH  J,B?-B1,I?  AS  FOLLOWS 

DURING  RUN  NO.  »**,  AH  WAS  IN  PLAY  FOR  •*.***  MINUTES, IZ=  * 

CALL  SCRIBE  (AHX  <*  ) , AHY  (» ) , ARTYX(»  ) , ARTYYt  • I , X .MAX , Y . MAX , K ,L , J , 1 1 , PK.  X (»  ) , PK  . Y « • ) 
,PKSI*> ,NO.RUNS,PE.R,PE.O,AIMX, AIMY.FLA5?) 

ALWAYS 
LET  J«J*1 

IF  J LE  NO.RUNS, SCHEOULE  A SCENARIO  NOW 

LET  1=0  LET  K*0  LET  L=0  LET  11=0  LET  FLAG=0 

LET  FLAG  1*0  LET  FLAG?*0 

FOR  EACH  AH  IN  POSITION, DO 

REMOVE  THE  AH  FPOM  POSITION 

DESTROY  THIS  AH 

LOOP 

UNTIL  N.EV.SCI.THOFAT>=0,D0 
LET  THREAT=F .EV.S(I.TPREAT) 

CANCEL  THREAT 
DESTROY  THREAT 
LOOP 

UNTIL  N.EV.S<I.PREC.FIRE)»0,00 
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LET  PREC.FIRE=F.EV.S< I. PREC.FIRE) 

CANCEL  PREC.FIRE 
OESTROT  PREC.FIRE 
LOOP 

UNTIL  N.EV.S<I.LOS)=0,DO 
LET  LOS=F.EV.S(I.LOS> 

CANCEL  LOS 
OESTROY  LOS 
LOOP 

RELEASE  AHX.AHf,  ART VX, A RTYY , PK. X, PK . Y , PKS 

RESERVE  ARTYX (*)  , ANO  ARTYYI*)  AS  999, AND  PK  .X  (*  ) , PK.  Y (*  > , ANO  PKSC)  AS  999 
SCHEDULE  A THREAT  IN  UN IFORH.F C 0 . 0 , T 1 ,5 > MINUTES 
IF  OPTION  EO  1, 

SCHEDULE  AN  LOS  IN  1./6.  MINUTES 

ALWAYS 

RFTURN 

ELSE 

PRINT  1 LINE  WITH  STAT  . 1 ,ST AT . ? AS  FOLLOWS 
AVG  P-K=  ».**♦♦>  SO=  *.«•#•* 

STOP 

END 

ENO  OF  RECORO  CARO 


USER’S  DATA  CARDS  ENTFREO  HERE 


6/7/S/9  ENO  OF  JOB  CARO 
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